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I. INTRODUCTION 



The Finite Element (FE) method is a proven tool for modeling structural dynamic 
systems. As the complexity of the system increases the FE model may not accurately 
reflect the dynamic behavior of the system To determine the extent to which the FE 
model accurately describes the physical system, a comparison of the dynamic response or 
modal parameters of the system as predicted by the FE model and the response or modal 
parameters of the physical system as determined by measurements of the dynamic 
response of the system can be made. Such a conparison can easily point out the 
differences in the dynamic behavior of the FE model and the physical system but feil to 
provide the necessary corrections to the FE model that will provide a more accurate FE 
model prediction of the dynamic behavior of the physical system 

Structural system identification refers to procedures used to identify finite element 
modeling errors using dynamic test data. Localization is the process of identifying those 
degrees of freedom (DOF) of the FE model whose impedance differs from that of the 
physical system We refer to this set of DOF as the error set. Identification is the process 
of finding matrices AAT, AM, and AC that are corrections to the FE model stiflftiess, mass, 
and damping matrices. When corrections are installed, the frequency response of the 
corrected FE model better predicts the frequency response of the experimental system, 
and hence the modal parameters of the corrected FE model better predicts those of the 
experimental system 

Structural system identification can be conducted using either modal or frequency 
domain methods. This thesis investigates a frequency domain method based on the 
structural synthesis transformation (SST) as outlined in Reference (1). When measiured 
dynamic response data is available for each DOF of each coordinate of the FE model the 
identification is termed spatially complete and the SST can be used to exactly determine 
which elements of the FE model are in error. Additionally, the SST can be used to 
conqjute correction matrices AM, AAT, and AC that can be used to correct the mass, 
strffiiess, and damping matrices of the FE model. The dynamic response of the corrected 
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FE model exactly matches that of the physical system. In the more frequent case that 
dynamic response data is not available at every DOF of the FE model, the identification is 
termed spatially incomplete. In this case the SST provides a frequency dependent 
solution. 

SST-based structural system identification uses the frequency response fimction 
(FRF) of the structure under consideration. The FRF is a quantitative description of the 
dynamic behavior of the structure. To obtain the FRF, known harmonic excitation forces 
are apphed and the resulting harmonic response of the structure measured. The ratio of 
excitation forces to response at a coordinate evaluated at each frequency in a specified 
bandwidth defines the FRF. 
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n. THEORY 



A. IMPEDANCE DESCRIPTION 

The impedance model of a given physical system can be defined by the relationship 
of structural displacement with respect to an applied force. 




(2.1a) 



The harmonic force and response vectors are denoted by "f and "x" respectively. 
These vectors and the impedance matrix, Z, are in general conplex-valued and fi'equency 
dependent. Subscripts "i" and "c" denote non-error and error DOF respectively. The 
superscr^t "a" denotes quantities calculated firom a FE (analytic) model. If the values 
were obtained firom experimental test data, the superscr^t would be "x". Thus, for the 
experimental model the inpedance relationship would be given by; 




(2.1b) 



In pratice, the elements of the experimental in^edance matrix of Equation (2.1b) 
are unmeasurable quanties. To see this we e)q)and row j of Equation (2.1b) to obtain 



f> 



j = W +^/2^2+- 



■+2jnX„ 



( 2 . 2 ) 



To measure an element z', of T, requires that we inpose a unit displacement at 

coordinate xi while physically holding all other other coordinates at zero displacements. 
This is not physically possible. Ass uming , for pmpose of definition, the availability of the 
experimental impedance matrix, the quantitative difference between the analytical and 
experimental systems, as a function of fi’equency, is described by the impedance error 
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matrix. It is defined by the difference between the analytical and experimental impedance 
matrices. The error impedance matrix relationship is defined as: 



o 


1 

o 




~ya 


^ic 




' 7 ^ 


^ic 


0 


AZ(Q)_ 




7^ 


N 




Z"" 


Z"" 



B. STRUCTURAL SYNTHESIS TRANSFORMATION 

Since the experimental impedance matrix, 2*, is in general unavailable, fi-equency 
domain structural synthesis is used to identify the impedance error matrix using FRF data 
exclusively. A structural synthesis transformation is constructed firom AZ of Equation (2.3) 
which encompasses the FE model errors. This transformation is applied to the finite 
element model to produce an experimental system FRF. 

The FRF relates structural response to applied excitation. Given a FE model with 
impedance matrix, Z“, the FRF, }f, is the matrix inverse of the impedance matrix Z" of 
equation (2.1a). For a static system (Q=0) the FRF is simply the flexibility matrix (inverse 
stifl&iess matrix). Using the notation of Equation (2.1) we may partition If as follows: 






K 


K 


K 




(2.4) 



Generally, the "c" response coordinates experience applied forces due to both error 
impedances and externally applied forces, whereas "i" response coordinates experience 
only externally applied forces, such that, 

f,=/r+ff (2.5a) 



and 






(2.5b) 



4 



Expanding Equation (2.4) and substituting the relations of Equation (2.5) yeilds 



the following relationships: 

I, = Kfr + H‘jr + (2.6a) 

= Kfr + KfT + (2-6b) 



If we include a copy of Equation (2.6b), in expanded matrix notation. Equation 
(2.6) reflects the three harmonic excitation terms to be considered, i.e.. 







Tja Tjo Tja 

^ii ^ic ^ic 




'/r' 




. == 


rja Tja tjo 

^ci ^cc ^cc 




/r > 


/c. 




Tja Tja j^a 

_^ci ^cc 




> 



Response coordinates "c" and "i", which are due to external forces, will hereafter 
be referred to as "e” coordinates, denoting their dependence on external force excitation. 
Consequently, the three excitation forces are condensed into two under the identity: 

{/.} = [[/"] [r*]]" (2.8a) 



{/.} = {X“} 



(2.8b) 



and Equation (2.7) reduces to 




(2.9) 



Equation (2.3) shows that the impedance error is defined as the difference between 
the analytic and eiqjerimental impedance models. Hence, a transformation is required 
which uses the FRF relationship of Equation (2.9) to generate a similar relationship for the 
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experimental system. The impedance error AZ provides the basis by which this 
transformation is developed. 



The impedance error matrix, AZ(Q,), is a function of frequency and satisfies: 



{/.} = -[AZ(Q)]{a:4 



( 2 . 10 ) 



where 



[az(q)] = [[A/s:] + yo[AC] - q^[am]] 



( 2 . 11 ) 



for Q the forcing frequency, j = and AK, AC, and AM stiflSiess, damping and mass 
error matrices comparable to those of the finite element formulation. The minus sign in 
Equation (2.10) reflects that the reaction forces imposed by inqjedance errors on the 
baseline model are being considered. Substituting the relationship 



f/e' 

u. 



7 0 

0 -AZ. 



Ife' 



( 2 . 12 ) 



into Equation (2.9) yields: 



TjO ttQ 

^ee ^ec 

Tja Tja 

^ce ^cc 



7 0 

0 -AZ 



7; 



(2.13) 



simplifying we get: 



K -Kaz 
K -CAZ 



7" 



(2.14) 



Expanding Equation (2.14) into two equations and using to denote a 
synthesized modified response results in 



x ; = HU - 



(2.15a) 



xi=hu-h:,azx: 



(2.15b) 
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Rearranging Equation (2.15b) produces 

\I * Hlsz\x: = HU 


(2.16a) 


Using the property of the frequency response function, 

= Kf, 


(2.16b) 




(2.16c) 


x] =[/ + //;aZ]'’x. 


(2.16d) 


Introducing Equation (2.16) into Equation (2.15b) results in 

xl = HU, - ^iA2[/ + //J AZ]'‘ X, 


(2.17a) 


= HU, - H‘„/^ I + //JAZ]"' HU, 


(2.17b) 


Once again we recall the property of a FRF, 

= H-u, 


(2.18a) 


Combining Equations (2.17b) and (2.18a) yields 

Hi =H-„- //;az(/+ 


(2.18b) 


Noting that 

[/+//iAZ)''=[(4Z-‘+K)AZ]'' 


(2.19a) 


and applying the matrix property 

«»][»])■'=[*]■ VI" 


(2.19b) 


we get that 

Ke=H°^~ + kVk 


(2.20a) 
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Replacing the superscript which denotes the structures's synthesized coupled 
response, with the superscript "x" to indicate the test system response we arrive at 



Equation (2.20b) is the structural synthesis transformation equation. When the 
experimental system FRF, If , is available the SST can be used to identify a frequency 
dependent impedance error matrix [AZ(Q)]. Additionally, using Equation (2.9), [AZ(t^)] 
can be decomposed into constitutent stiffaess, mass, and damping errors. 

C. FREQUENCY DOMAIN LOCALIZATION 
We may rewrite Equation (2.20b) as 




(2.20b) 



In full matrix notation we have 




(2.20c) 




( 2 . 21 ) 



where 




(2.22a) 



and 






(2.22b) 



We define the localization matrix L as 




(2.23) 



using the expression of Equation (2.22a) in Equation (2.23) we can rewrite L as 



(2.24a) 



Expanding the "e" coordinate set into error and non error coordinates we get 



8 



(2.24b) 




Noting the the frequency response matrix is the inverse of the impedance matrix, i.e., 





Z""’ 




K' 




7 


o' 


\_^ci 


1 

N 


K 


Kc. 




_o 


7 



(2.25a) 



all mixed product coordinates of Equation (2.24b) must be zero and L simplifies to 



L = 



"0 0 ■ 
0 /)■’ 



@Q=Qi 



(2.25b) 



D. ERROR IMPEDANCE 

We can solve Equation (2.20b) for the impedance error [AZ]. From Equation 
(2.25b) terms of Equation (2.20b) associated with non error coordinates may be assumed 
to be zero. We get the following form of the inpedance error matrix 

[A2]=([#;]-[/fj])“ @n=£ii (2.26a) 

where 

= . (2.26b) 

Let H={Qi, Q 2 , ... , On} be a set of frequencies where Qi < Q 2 <...<On-i<On. If 
for each i= 2,3, ..., n-1, we apply Equation (2.10) at each of the frequencies Qj.i, Qj , and 
Qi+i and assemble the resulting three equations into a system of three equations in three 
unknowns, we get the matrix equation 



■aZ.(£2,_,)' 




1 

D 




'AK/ 


AZ,(n,) 


= 


I -Cl]l 


jClJ 


AM, 


.AZ.{n„,) 




J 







Equation (2.27) can be used to decompose the frequency dependent impedance 
error into constitutent stiffiiess, mass, and damping error matrices, AK, AM, and AC at the 
frequencies Qj, i=2,3,...4i-l. These constitutents matrices are in general frequency 
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dependent. Denoting the solution of Equation (2.27) by 



AC(n,) 



for i=2,3,...,n-l we 



obtain for each Qi i=2,3,...,n-l, error stifl&iess, mass, and damping matrices AKc(Qi), 
AMc(f2i), and ACc(Q) such that 



AZ,(n,) = - n? am.(£1,) + JdAC^n,) 



(2.28) 



The matrices AKc(Qi), AMc(f2i), and ACc(f2i) can be used to numerically correct 
the stifl&iess, mass, and damping matrices of the FE model at the frequencies Qj 
i=2,3,...,n-l in that the FRF of the corrected FE model at the frequencies Q; approximates 
the experimental system FRF at Qi. To express this symbolically, if AT”, Kf, and C are the 
stiffness, mass, and dampit^ matrices of the FE model and If is the FRF matrix of the 
experimental system at Qi 

» (k- + AK,{n,) - nf (at + AA4{n,)) + ;n,(c” + AQ(n,)))‘' (2.29) 



where the "c" subscripted matrices are added at the corresponding error set coordinates of 
the "a" superscripted matrices. 

For a general set of frequencies, S={Qi, Q2, ... , Qn}, we can form the system of n 
equations in three unknowns given by: 



‘AZ.(n,)‘ 




'/ -Q^/ ;Q,/' 


AZ.(£2,) 


= 


/ -Q^7 JQ,I 


.AZ.(n.)_ 




J -all yQ„/J 



AM, 

AC, 



(2.30) 





'AK,' 




■ak.(=)' 


The solution 


AM, 


= 


AM,(=) 




ACc. 




_AC,(H)_ 



of Equation (2.30) represents error stiffiiess. 



mass, and damping matrices which best corrects the FE model in a least squares sense. 
Equations (2.27) and (2.30) are fundamental to all that follows. 
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m. SPATIALLY COMPLETE STRUCTURAL IDENTIFICATION 



To illustrate the principles of SST based frequency domain structural identification 
we will use the FE model of a free-free beam. To simulate the experimental system we 
impose a 25% addition to the mass and strfiBiess of elements 3 and 4 of a 10 element FE 
model. Figure 3-1 shows the finite element model of the beam and the spatially complete 
ejqjerimental system that results from imposing the mass and stifiEhess additions at element 
3 and 4 of the FE model. Table 3-1 shows the system frequencies of the analytic and 
ejq>erimentai systems. 



7^: 7 


7 




7 


7 


7 


7 


7 


7 


“s 7 




1 ^ 




1 ^ 


1 ^ 


1 ^ 




1 ^ 




lJ 



EXPERIMEKTAL SYSTEM 
MASS OF ELEMENT 3 * 1.954X10^ UBM 
MASS OF ELEMENT 4 » 1.954X1(r-6 LBM 
STIFFNESS OF ELEMENT 3 * 417.88 LBF/)N 
STIFFNESS OF ELEMENT 4 » 417.88 LBF/IN 






/N 



/N 



/In 



/N 



/K /K /N /K /K /N 









1 ^ 


1 ^ 




1 ^ 






1 ^ 



ANALYTIC SYSTEM 
NUMBER OF ELEMENTS » 10 
LENGTH *60.625 
EJ* 74500 LBF4N*2 
SECTIONAL AREA * 0.8359 IN*2 
WEIGHT DENSITY * 0.05 LBF/IN*3 
ELEMENTAL MASS * 1.563X10^6 LBM 
ELEMENTAL STIFFNESS * 334.3 LBF/IN 




Figure 3- 1 Spatially complete Analytic and E^erimental Systems. 



11 



MODE 


ANALYTIC (Hz) 


EXPERIMENTAL 

(Hz) 


1 


25.41 


24.42 


2 


70.06 


68.41 


3 


137.4 


130.2 


4 


227.5 


216.3 


5 


340.8 


329.4 


6 


478.0 


456.5 


7 


639.9 


612.0 


8 


826.5 


802.6 


9 


1026 


969 


10 


1363 


1323 


11 


1641 


1573 


12 


1981 


1920 


13 


2381 


2283 


14 


2850 


2754 


15 


3397 


3269 


16 


4028 


3912 


17 


4720 


4605 


18 


5377 


5156 


19 


6795 


6791 


20 


6808 


6801 



Table 3-1 System Frequencies of spatially complete Analytic and Experimental systems. 

We will denote by M “ , K ° , and C" the mass stififiiess and dancing matrices of 
the FE model and hyM^, K"‘ , and the mass stififiiess and damping matrices of the 
experimental system. The impedance matrices of the analytic and simulated experimental 
systems are given by: 
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(3.1a) 



Z‘‘{a) = K“ +jQ.C“ 

Z^(Q) = K^ + jQC" - (3.1b) 

The FRF matrices of the analytic and simulated experimental systems are given by: 

H“{n) = z^ny^ (3.2a) 

W{Q) = Z^(Q)’’ (3.2b) 

Figure 3-2 shows a comparison of the driving point FRF at DOF 1 of the analytic 
and experimental system. 
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Omega (Hz) 



Figure 3- 2 Analytic FRF vs simulated Experimental FRF. 
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• For a given frequency Qo, using Equation (2.24a), we can form the localization 
matrix L. Plotting the diagonal elements of L versus the associated DOF we obtain the 
plot shown in Figure 3.3. For each frequency, Q, in a given frequency range, we can 
compute the diagonal of L at . Assembling the diagonals over the frequency range of 
our system into a rectangular matrix and performing a MATLAB mesh plot of the 
resulting rectangular matrix we obtain the surface plot shown in figure 3.4. Figure 3.4 
shows the frequency dependency of the localization matrix diagonals. As our system is 
spatially conq)lete Equation (2.25) forces all non error coordinates to be zero. From 
Figure 3-3 we can determine that the locations of the nonzero diagonal values are DOF 9, 
10, 11, and 12. We denote by the set of DOF for which L(i,i) is non zero. For our 

beam system = {9,10,11,12} . Figure 3.5 shows the frequency dependence of typical 
error and non error set diagonal elements of the localization matrix L. 

Using the set, Cetr, which results from the localization, we can perform a 
partitioning of the FRF matrix as described by Equation(2.4). We can now apply 
Equation (2.26) to compute AZ as a matrix fimction of frequency over the frequency 
range of our system. We then use Equation (2.27) to decompose AZ into its constituent 
components AK, AM, and AC. We get exact solutions of error stiflhess, mass, and 
damping as shown in figures 3-7, 3-8 and 3-9. The MATLAB Routine SST.M of 
Appendix A can be used to acconplish the above steps. 

For each Q in the frequency range of our system we can form the sum: 



as the corrected FRF of the analytic model. Hcon- is the FRF of the model that results from 
installing the corrections as identified by the Equation (2.26a). Figure 3-10 shows a 
comparison plot of the FRF of the corrected model and the FRF of the experimental 
system. Figure 3-10 clearly shows the exactness of the SST solution in the case of a 



z^(n) = zi{n)+e^(0) 



(3.3) 



We refer to 




(3.4) 
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spatially complete system. The experimental and corrected model FRPs are identical to 
within plot resolution. 

At each frequency, Q, in the frequency range of our system we can form the sum 
Z^(Q) = (Q) + A^(Q) + yQAC(Q) - Q^AA/(Q) (3.5) 

We refer to: 

(3-6) 

as the constituent corrected FRF. Figure 3-11 is a comparison plot of the constkutent 
corrected FRF and experimental for our system. The offset between the two plots is 
exactly equal to the sampling frequency AQ. 
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Figure 3- 3 Spatially complete localization matrix diagonal at Q=352 Hz. 
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Figure 3- 4 Frequency dependence of spatially complete localization matrix diagonals. 
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Figure 3- 5 Frequency dependence of spatially complete localization matrix error set 
DOF. Units are IbFin (LoglO of). 
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Figure 3- 6 Spatially complete Analytic Impedance vs Experimental impedance. 
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Figure 3- 7 Computed stifi&iess vs true stiffiiess for spatially complete beam. Plots are 
identical within plot resolution. 
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Figure 3- 8 Computed mass vs true mass for spatially complete beam. Plots are identical 
within plot resolution. 
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Figure 3- 9 Computed damping vs true damping for spatially complete beam. Plots are 
identical within plot resolution. 
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Figure 3- 10 Spatially complete Experimental FRF vs AZ corrected FRF. Plots are 
identical within plot resolution. 
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Figure 3-11 Spatially complete experimental FRF vs AK, AM, and AC corrected FRF. 
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IV. SPATIALLY INCOMPLETE STRUCTURAL 
IDENTIFICATION 



A. GENERAL DESCRIPTION 

To illustrate SST based frequency domain structural identification when the 
physical system imder consideration is spatially incomplete we will ^ain use the FE model 
of a free-free beam. We simulate the experimental system by imposing a 25% addition to 
the mass and stififtiess of elements 3 and 4 of a 10 element FE model. Figure 4-1 shows the 
FE modeled beam and the spatially incomplete system that results if FRF data is available 
only at the displacement DOF of the FE system. 
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ANALYTIC SYSTEM 
NUMBER OF ELEMENTS = 10 
LENGTH =00.625 
B = 74500 LBF-IN*2 
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Figure 4- 1 Analytic and experimental spatially incomplete systems. 
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To obtain a simulated FRF for our spatially incomplete beam we denote by 
M° , K° , and C° the mass stifi&iess and damping matrices of the FE model and by M"' , 
AT"' , and the mass stiflBiess and damping matrices of the experimental system. , 

K "" , and are obtained by imposing 15% mass and stififiiess and a 15% mass additions 

to the elemental matrices of element 5 and 6 respectively of the FE model. The impedance 
matrix of the simulated spatially complete beam is given by 

(Q) = + jCia - (4- 1 ) 

The FRF matrix of the simulated spatially complete system is given by 

//"(Q) = Z"(Q)'* . (4-2) 

We introduce the terminology Analysis set and Omitted set where the Analysis set 
(A-set) is that set of DOF for which experimental FRF data is available and the Omitted 
set (0-set) is that set of DOF of the experimental system for which experimental FRF data 
is unavailable. For our simulated experimental system the A-set consist of the odd 
numbered translational DOF and the 0-set consists of the even numbered rotational DOF, 
i.e., 

^-^er = {l,3,5,...,21} (4.3a) 

0-5c/ = (2,4,6, ...,22} (4.3b) 

We obtain the simulated FRF of our spatially incomplete beam by physically 
extracting the rows and columns of the simulated spatially complete matrix. If, for which 
FRF data would be available. In our system these are the rotational DOF and all even 
numbered rows and columns are omitted from the simulated spatially complete FRF to 

obtain a spatially incomplete FRF which we will denote by H "‘ . 

For a fixed Qo, is a square matrix of size (length(yf-sc/)) by (length(.4- 

set)). For our FE model as currently defined. If, is of size (number of DOF) by (number 

of DOF) and is, as is true of most real world cases, of larger size than . In order to 
employ the structural synthesis transformation we need to reduce the size of If to that of 
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. To this end we will consider two reduction methods, FRF matrix extraction [Ref. 1], 
and the Improved Reduced System as given in [Ref. 2]. 



B. EXTRACTION REDUCTION METHOD 



In order to reduce ff by the extraction method we simply extract from the full 
order H those rows and column which correspond to A-set coordinates. Partitioning the 
impedance and full FRF matrices of our analytical system 




(4.4a) 



H = 



Tja ttQ 

^aa ^ao 



Tjo j^a 



(4.4b) 



and using the identity 



-7a 


1 

N 




HI 




7 


O' 


_^oa 


^00 _ 


HI 


Hl_ 




_0 


/_ 



(4.5) 



we obtain the relationship 






(4.6) 



We denote the reduced Analytic FRF of Equation (4.6) by H° . 



C. 0-SET SYSTEM 

Equation (4.6) relates the extracted reduced order FRF of the analytic model to the 
impedance of the full order order model. Taking advantage of the identity 

[Z.„]"'=(det[Z„])‘'arf;[Z„] (4.7) 

and replacing Zoo by Koo+j^2Coo-J^Moo (where det(*) and adj(*) represent the 
determinant and adjount respectively), we see that an element //^(Q) is large for those 

frequencies Qo where Qo is an eigenvalue of the 0-set system, the 0-set system being that 
FE model having stiffiiess, mass and damping matrices Koo, Moo, and Coo respectively. 
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D. IMPROVED REDUCTION SYSTEM 



To use the improved reduction method (IRS) we first partition using the A and 
O sets, then adjust Equation (2. 1) to reflect this new coordinate system obtaining 




Expanding Equation (4.8) into two equations yields 

X = Zlx, + 


(4.9a) 


f X ^2° X 

Jo ^oa'^a~ ^oo'^'o 


(4.9b) 



0-set coordinates are not associated with FRF data measurement locations on the 
physical structure, therefore the forcing fimction at 0-set coordinates can be set to zero. 
Making these substitution into Equation (4.9b) and solving for the generalized structural 
response coordinates leads to: 






° X 
oa’^a 



(4.10a) 





(4.10b) 



Substituting these results into Equation (4.8) yields, 





^oo 



rya rya 

^oo ^oa 




(4.11) 



hence 

{/.}=[z«-z»z»z»lk} (‘*■12) 



When 0=0 the Equation (4.10a) yields the static reduction relationship between 
omitted and retained coordinates and is given by 

(4.13) 
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The IRS relationship is given by 






(4.14) 



where 




(4.15) 



and Kstat and Msut are the statically reduced [Ref. 2 , 3] stiffiiess and mass matrices. 

Unlike the spatially complete case where we only had to consider two system (the 
analytic and the experimental systems), in the case of a spatially incomplete system there 
are actually five systems with which we must concern ourselves [Ref 4]: the analytic 
system, the experimental system, the reduced analytic system that results from conducting 
dynamic reduction on the mass and stiffiiess matrices of the analytic system, and the 
omitted systems of both the analytic and the experimental system. Table 4-1 shows the 
frequencies of each the first four of these systems. 

Figure 4-2 show^ a comparision of the analytic and experimental FRF of our 
spatially incomplete beam. We see that by using IRS reduction of the analytic system those 
modes above approximately 1000 Hz i.e., those modes associated with the reduced out 
rotations, are not present. Figure 4-3 shows a similar comparision where we have used 
extraction reduction and the higher modes are present. In all that follows we shall use IRS 
reduction of our analytical systems. Figure 4-4 shows the localization matrix diagonal at 
Q= 196.1 Hz while Figure 4-5 shows the frequency dependency of the localization 
diagonals over the frequency range our our spatially incomplete beam. For a spatially 
incomplete system the determination of the locations of the error coordinates is not a 
clearly defined task. We shall not discuss the problem of actually determining the exact 
error set in the spatially incomplete case and will use our knowledge of the true location of 
the error coordinates to aid in our localization. As we shall make use of the concept of the 
size of the error set again, we will simply note that using a reduction method like IRS 
causes errors to be 'smeared' in the reduced analytic model. Table 4-2 shows the A-set and 
0-set of our simulated beam along with the locations of the true errors as well as the 
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computed C-set at Q=196 Hz. The computed C-set is the set of all DOF having a diagonal 
entry whose absolute value exceeds a given tolerance. 



Mode 


Anal (Hz) 


Ejq) (Hz) 


Reduced 

(Hz) 


0 set (Hz) 


1 


25.41 


25.83 


25.41 


1244.28 


2 


70.07 


70.85 


70.07 


1287.94 


3 


137.45 


137.93 


137.45 


1417.14 


4 


227.55 


227.37 


227.55 


1629.59 


5 


340.84 


341.14 


340.84 


1928.6 


6 


478.0 


477.97 


478.0 


2327.83 


7 


639.9 


639.92 


6403 


2853.66 


8 


826.5 


826.24 


829.7 


3541 


9 


1026.8 


1025.99 


1029.74 


4399.32 


10 


1363.62 


1364.51 




5280.86 


11 


1640.62 


1640.71 




5701.99 


12 


1981.19 


1984.61 






13 


2381.4 


2380.32 






14 


2850.01 


2847.38 






15 


3397.39 


3397.74 






26 


4027.82 


4026.28 






17 


4720.15 


4738.64 






18 


5376.83 


5380.09 






19 


6794.91 


6787.09 






20 


6807.52 


6804.12 







Table 4-1 Analytic, Experimental, Reduced, and Omitted System Frequencies of Spatially 
Incomplete Beam- 
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A Set 


1 


3 


5 


7 


9 


11 


13 


15 


17 


19 


21 


0 Set 


2 


4 


6 


8 


10 


12 


14 


16 


18 


20 


22 


True C Set 


5 


6 


7 


8 


9 


10 












Computed C 
Set 


1 


3 


5 


7 


9 


11 


13 











Table 4-2 Analytic set, Omitted set, Computed C set, and True C set DOF for spatially 
incomplete beam. 
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Figure 4- 2 Anal 5 ^ic FRF vs experimental FRF for spatially incomplete beam using IRS 
reduction. 
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Figure 4- 3 Analytic FRF vs experimental FRF for spatially incomplete beam using 
extraction reduction. 
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Figure 4- 4 Localization matrix diagonal at Q=196. 1 Hz for spatially incomplete beam. 
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Figure 4- 5 Frequency dependence of localization matrix diagonals for spatially 
incomplete beam. 
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Figure 4- 6 Frequency dependence of error and non error set localization matrix 
diagonals for spatially incomplete beam. Units are Ibfi'in (Log 10 oQ. 
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Figure 4- 7 Analytic vs experimental impedance for spatially incomplete beam. 
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Figure 4- 8 Computed Stifi&iess vs True StifBiess for spatially incomplete beam. 
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Figure 4- 9 Computed Mass vs True Mass for spatially incomplete beam. 
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Figure 4- 10 Computed Damping vs True Damping for spatially incomplete beam. 
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Figure 4- 11 AZ Corrected FRF vs experimental FRF for spatially incomplete beam. Plots 
are identical to within plot resolution. 
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Figure 4- 12 AK/AM/AC Corrected FRF vs experimental FRF for spatially incomplete 
beam. The offset of the two plots is equal to the sampling frequency AQ. 
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V. SINGLE MODE SOLUTIONS 



A. SINGLE MODE MATRIX SOLUTIONS 

In the case of a spatially complete system we have seen that the SST yielded 
frequency independent error matrices AK, AM, and AC that could be used to correct the 
stiffiiess, mass, and damping matices of the FE model in such a manner that the FRF of the 
corrected FE model was exactly equal to that of the experimental system. In the case of a 
spatially incomplete system the constituent solutions AK, AM, and AC given in chapters 
III and IV were in general frequency dependent solutions which serve only as corrections 
to the reduced FE model. Ideally we seek frequency independent solutions which are 
corrections to the full FE model. For now we shall only deal with the simpler problem of 
trying to find frequency independent solutions which serve as corrections to the reduced 
FE model. 



We shall employ Equation (2.30). For a given experimental system mode, (On, 
consider a frequency bandwidth [Qi, Qu] such that tOne[Qi, Qu]. Let S={f^i Q 2 ,...,f^m} be 
a frequency sampling of the bandwidth [Qi, Qu], i.e., Qi<Qi<Qufor each QicH. For each 
QiGH we can form the error impedance matrix AZ(/2i). We apply Equation (2.30) to the 
partition yeilding the following system of m equations in three unknowns: 









= 


_AZ.(£2.)_ 





-Q(7 yQ,7 
-Qf7 ;q,7 

j^J 






AA7;» 
ACf- 



(5.1) 



We will demonstrate by example that for properly chosen bandwidths [Qi, Qu], the 
solution AAT®", AM®”, AC®" of Equation (5.1) approximately corrects the Analytic FRF 
over the frequency bandwidth [Qi,Qu], i.e., if ff(jO) is the value of the FRF matrix of the 
experimental system at a frequency Q, 77‘'(Q)the value of the FE reduced analytic FRF 
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matrix at Q, and H°„^{Ci) the value of the single mode corrected reduced Analytic FRF 
matrix at Q for an experimental system mode cOn where Qe [f2i,f2u] then 



H‘(0) - /C,(£2)| S |W'(£2) - H"(n)| 



(5.2) 



We will refer to the solution AAT®" , AA/"" , AC“" as a single mode solution at cOn- 

In the case of a spatially complete system, such as the spatially complete beam discussed in 
chapter III, all single mode solutions are found to be identical to the unique frequency 
independent solution that was obtained in chapter III, hence the corrections are vaUd 
throughout the frequency range of the experimental system and the left hand side of 
Equation 5.2 is zero. For a spatially incomplete system this is not the case. 

In support of Equation (5.2) we chose a 1 Hz frequency bandwidth centered on 
mode 1 (coi=25.2I Hz) of our spatially incomplete beam as defined in Chapter III. For our 
frequency sampling, E, we will use a sampling frequency of .5 Hz which results in 3 points 
which are equally spaced over the interval, E={24.7178 Hz, 25.2178 Hz, 25.7178 Hz}. 
Solving Equation (5.1) we get 



= 



-3.91 

-13.60 

-21.70 



-13.60 

170.5 

-187.5 



21.70 

-187.5 

194.6 



(5.3a) 



AAC' = 



- 0.0001 

0.0006 

-0.0005 



0.0006 

-0.0033 

-0.0026 



-0.0005 

0.0026 

-0.0016 



Acr = 



0 + 0.0270y 
0-0.229iy 



0 + 0.2277y 



0-0.229iy 
0 + 1.5078y 
0-1.4627/ 



0 + 0.2277/ 
0-1.4627/ 
0 + 1.3781/ 



(5.3b) 



(5.3c) 



Figure 5-1 shows a comparison of experimental, tmcorrected and mode 1 
corrected FRFs for our spatially incomplete beam using a 1 Hz frequency bandwidth 
centered on Mode 1 and a frequency sampling consisting of 3 frequencies equally spaced 
over the bandwidth. Figure 5-2 shows the results of including mode 1 and mode 2 in the 
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frequency bandwidth [Qi, Ou]- As Figure 5-2 shows Equation 5.1 is not as accurate when 
multiple frequencies are included in the bandwidth. Figure 5-3 is a comparison plot of the 
experimental and analytic FRF versus the single mode matrix solution corrected FRFs over 
25, 10, and 1 Hz bandwidths using a 3 point frequency sampling. Figure 5-3 shows that 
the precedure is sensitive to the size of the bandwidth over which the solution is 
computed, i.e., better accuracy is achieved with smaller bandwidths. Figure 5-4 is a 
comparison plot of the experimental and uncorrected analytic FRFs versus the single mode 
matrix solution corrected FRFs computed over a 1 Hz bandwidth using frequency 
samplings of 200, 50, 10, and 3 points equally spaced over the bandwidth. Figure 5-4 
shows that the procedure is friirly insensitive to sampling frequency. 

Figure 5-5 is a comparison plot of the experimental and analytic FRFs versus the 
single mode solution corrected FRF for a 1 Hz bandwidth using a 3 point frequency 
sampling in the case of a spatially complete beam. As Figure 5-5 shows the procedure is 
exact for spatially complete systems. 
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Figure 5- 1 Experimental and imcorrected Analytic FRFs vs single mode solution at 
mode 1 corrected FRF using a 3 point frequency sampling of a 1 Hz bandwidth. 
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Figure 5- 2 Experimental and uncorrected aneilytic FRFs vs single mode corrected FRF 
using a 15 point frequency sampling of a bandwidth that includes modes 1 and 2. 
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Figure 5- 3 Spatially incomplete ejqperimental FRF vs single mode matrix solutions at 
mode 1 using 3 point frequency samplings of 25, 10, and 1 Hz bandwidths. 
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Figure 5- 4 Spatially incomplete experimental FRF vs single mode matrix solutions at 
mode 1 using 3, 10, 50, and 200 point frequency samplings of a 1 Hz bandwidth. 
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Figure 5- 5 Spatially complete experimental FRF vs single mode matrix solutions at 
mode 1 using a 3 point frequency samplings of a 1 Hz bandwidth. 



51 



B. SINGLE MODE INTEGRAL SOLUTIONS 



In what follows we wish to employ the integral formulas of reference (5) to obtain 
a results similar to that of Equation (5.2). To accomplish this we shall need to recast the 
in^edance equations of Chapter II in terms of velocity. We start with the equation of 
motion of a general 2nd order linear system 



Kx + Cx + Me = / 


(5.4a) 


where 

x=Xe^^ 


(5.4b) 


II 


(5.4c) 


If we differentiate Equation (5.4b) we get that 

X = jQXe^ = jOx 


(5.5) 


hence 

X 


(5.6) 


jO 



differentiating Equation (5.4b) twice we get that 

X = jd^jOXe^) = jOx 



(5.7) 



substituting Equations (5.6 )and (5.7) into Equation (5.4a) we get 

K— + Cx + MjQx = / (5.8) 

jCi 



We can rewrite equation (5.8) as 



— + C + MJQ]x=f 

JO J 



(5.9) 



Following reference (5) we can take the Laplace transform of equation (5.9) and 
thus write the impedance of the linear system in terms of the complex Laplace parameter s 
as 
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m, 

■’ d 



(5-10) 



Z{s) = Ms + C + — 
s 



In reference (5) the impedance matrix of a general system is represented as an 
infinite Laurent series about the orgin in powers of the complex Laplace parameter s as 

Z{s) = = Um,^AAs’+-+As+A, + A, -+■■■■<■ A.,\) (5.11) 



After defining a truncated Laurent series 



Z{s) = A^s+Ao+A^- = Ms + C+ — 
s s 



(5.12) 



which approximates the impedance matrix, an error fimction E(s), and a cost J are defined 



as 



and 



E{s) = Z{s)-Z{s) 






(5.13) 



(5.14) 



where fV(sJ is a complex valued weighting fimction, the subscript E denotes the euclidian 
norm of a matrix and the integration is performed over a prescribed path P in the complex 
plane. By setting the partial derivatives of J with respect XoM ,K , and C to zero the 
authors of reference (5) obtained expressions for matrices Af, , andC which 
approximate the stifi&iess, mass and damping matrices about a mode of the system. The 
e5q)ression for M ,K , and C are as follows: 



M= ^ 



ab-c 



b |qZ;(/Q)|IT(/Q)|c/Q - c \—Zi{in)\W{iCl)\dCl 

Qi £1,^ 

(5- 15a) 

= -^Z^{i€t)\W{iCl)\dCl (5-15b) 
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where 



K = 



1 

ab-c^ 



O2 ^1 

c \QZi{iCL)\W{iCl)\iCl-a \—Zj(iCl)\W{iCl)\dO. 
. fi. 



(5-15C) 



^2 

a= |Q^|r(/Q)|i/Q 
*= 

«> 

c= ||rr(/Q)|c?Q 

«■ 

Z(7Q) = Zn{id) + /Z,(/Q) 



(5-15d) 

(5-15e) 

(5-15f) 

(5-15g) 



We will apply Equation (5.15) to AZ as defined by the SST and obtain 
matrices AAT , AM , and AC which will serve as correction matrices of the FRF of the 
system about a given mode of the system. 

As an example, let us use Equation (5.15) to compute the (1,1) element of the 
correction matrix AC at mode 1 (®i=25.21 Hz) for our spatially incomplete system. We 
take the weighting fimction to be W(iQ)=l/(iQ) , the path P is taken to be along the 
imaginary axis fi’om 155.30 rads/sec to 161.58 rads/sec. 

To compute this line integral we shall use the trapezoid rule with a three point 
fi’equency sampling 

E = {155.30/*a<5fe / sec,\5SA4rads / sec,161.58ra<a[y / sec} (5.16) 

of the straight line path P (sampling firequency of % rads/sec). We first compute the value 
of the constant a of Equation (5 . 1 5d). Using vector notation we have that 

Q(E) = [155.30rad:s / sec,158.44ra7iy / sec,\6l.5Srads / sec] (5.17a) 

and 

1F(H) = [0 - 0.0064 j sec/ rads Si - 0.0063 j sec/ rads Si — 0.0062 j sec/ rads^ (5. 1 7b) 
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(5.18) 



hence the integrand is the vector 

Q^(E)-|)!r(E)| = [155.30,158.44,161.58] 

Using the trapezoid rule to compute the integral, we have that 

<1 = 7’ n’(S)W3)|<a = |fl^|^ + 15&44 + l^) =995.56™*)^; (5.19) 

155.30 II L 

b=.00000158 sec^/rad^ and c=0.0397 are computed in a similar manner. 

To actually compute the matrix AC ,we use the SST to confute the matrices 
AZ(155.30 rads/sec), AZ(155.44 rads/sec), and AZ(161.58 rads/sec) for the points of the 
frequency sam pling of the straight line path P. Using Equation (2.26) the matrix 
AZ(1 55.30 rads/sec) is seen to be 



AZ(l55.44™%,) 



(0.64£-9) + 0.02y 
(-0.4E-8)-0.04y 
(0.3E-8) + 0.0iy 



(-0.04£-8)-0.04y 

(-0.4£-9)-0.04y 

(0.25E-7)-0.10y 



(0.03£-7) + 0.01y 
(0.25E-7)-0.10y 
(-0.24£-7) + 0.15y_ 
(5.20) 




The integral in Equation (5.15b) is computed on an element by element basis. For 
the (1,1) element of Equation (5.15b), if we collect the (1,1) elements of the matrices 
AZ(155.30 rads/sec), AZ(155.44 rads/sec), and AZ(161.58 rads/sec) and use the 
weighting vector given in Equation (5. 17b) the integrand of Equation (5.15b) is the vector 



AZ^*'‘>(E)|1F(E)| = [0.0064 0.040 -0.848][0.64 0.63 0.62](£:-9)^*-%jSe^^^^ 

(5.21a) 

AZr(H)|)^(H)| = [0.0041 0.2885 0.0036X£-9)^*-%jSe^^^ (5.21b) 



Computing the integral we have 



AC(U) = 0.0232£ - 6 



O'MU 02885 . 






in 
(5.22a) 

(5.22b) 



Performing the above procedures for the remaining elements of Equation (5.15b) we get 
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AC = 



(5.23a) 



0.0232 -0.0443 

-0.0443 0.1278 

0.0300 -0.0914 



0.0300 

-0.0914 

0.0645 




In a similar manner we have 



AM = 



- 0.0001 

- 0.0001 

0.0002 



- 0.0001 

0.0015 

- 0.0020 



0.0002 

- 0.0020 

0.0027 



Ibm 



AK = 



-6.0476 

4.5347 

3.6757 



4.5347 

51.2196 

-71.7061 



3.6757 

-71.7061 

85.4417 




(5.23b) 



(5.23c) 



Figure 5-6 is a comparison plot of experimental and uncorrected analytic FRFs 
versus single mode integral solutions at mode 1 using a frequency sampling of 3 points 
equally spaced over bandwidths of 25, 10, and 1 Hz. As with the matrix solutions, the 
integral solutions are sensitive to the length of the bandwidth over which the solutions are 
computed with smaller bandwidths yielding more accurate solutions. Figure 5-7 is a 
comparison plot of experimental and uncorrected analytic FRFs versus single mode 
integral solutions at mode 1 over a 1 Hz bandwidth using frequency samplings of 3,10,50, 
and 200 points from the bandwidth. As with the matrix solutions, the integral solutions are 
also insensitive to the sampling frequency. 



Figures 5-8, 5-9, and 5-10 are comparison plots of matrbc and integral solutions at 
mode 1 computed over 25, 10, and 1 Hz bandwidths respectively, using three point 
frequency samplings. These three plots show that for these frequency sampling, the matrix 
solutions are more accurate than the integral solutions. We must note that we have used 
the trapezoid rule for computation of the integrals. It is expected that better accuracy 
would be achieved from the integral solutions if a higher order integration method such as 
Simpson's rule were used. 
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Figure 5- 6 Spatially incomplete ejqperimental FRF vs single mode integral solutions at 
mode 1 using 3 point frequency samplings of 25, 10, and 1 Hz bandwidths. 
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Figure 5-7 Spatially incomplete experimental FRF vs single mode integral solutions at 
mode 1 using 3, 10, 50 , and 200 point frequency samplings of a 1 Hz bandwidth. 
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Figure 5- 8 Experimental FRF vs single mode integral and matrix solutions at mode 1 
using a 25 Hz bandwidth with a 3 point frequency sampling. 
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Figure 5- 9 Experimental FRF vs single mode integral and matrix solutions at mode 1 
using a 10 Hz bandwidth with a 3 point frequency sampling. 
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Figure 5- 10 Experimental FRF vs single mode integral and matrix solutions at mode 1 
using a 1 Hz bandwidth with a 3 point frequency sampling. 
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VI. MULTIPLE MODE SOLUTIONS 



A. MULTIPLE MODE MATRIX SOLUTIONS 

We wish to extend the results of Chapter V to multiple modes of the system under 
consideration. To accomplish this we shall simply extend Equation (5.1) over multiple 
frequency samplings. For a set of system modes {(Oi, C 02 v>®n}, let S(cOi)={Qii, 
Oi 2 ,.-MOik} be a frequency sampling of a bandwidth [Qu,Qiu] about ©i for i=l,2,...,n. We 
shall apply Equation (2.30) to the frequency sampling 

which is the set theoretic union of the sets H(©i) for i=l,2,...,n. If, for simplicy, we restrict 
ourselves to equally sized frequency samplings we obtain a set on nk equation in three 
unknowns 



■aZ.(£5„)' 




■/ -Q^,/ 


J^l\I 






AZ.(n,) 


— 


/ -Qj./ 


jW 


'AKf 

AMf 


(6.2) 






• • 


: 












J^nk^ _ 







One can solve Equation (6.2) for 



ak: 

am; 

ACf 



as we have done in Chapter V using the 



familiar MATLAB puedoinverse fimction but better results are obtained if we weight the 
equations as discussed in [Ref 6]. A good weighting is to assign the weight coi to those 
equations associated with the frequency sampling H(©i). Alternately we could assign the 
weight l/( 0 i to those equations associated with the frequency sampling 5((Di). The ©i 
weighting results in the solution being more accurate at the lower modes while the 1/©, 
weighting gives better accuracy at the higher modes. 



We note that the solution to Equation (6.2) approximately corrects the analytic 
FRF in the sense of Equation (5.1) when the matrices AZ(Q\^ as given by Equation (2.26), 



62 



are computed over a 'good' error set, Cea. In general, due to the smearing that occurs 
during reduction of the full analytic system to the reduced analytic system, the set Cerr is 
not the same set as the intersection of the spatially complete error set and the A set. It is in 
general a larger set than this intersection and the size is a function of the type of error, i.e., 
mass, damping, or stififiiess error and the locations of the errors. Figure 6-1 is a 
comparison plot of experimental and uncorrected analytic FRFs versus the multiple mode 
matrix solution corrected FRF with the solution having been computed over modes 1 and 
2 using a 1 Hz bandwidth about each mode with a 3 point frequency sampling over the 
bandwidth at each of the modes. Figure 6-2 is a comparison plot using a matrix solution 
computed over mode 1 through 4 again using a 1 Hz bandwidth about each mode with a 3 
point frequency sampling over the bandwidth at each of the modes. 
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H(9,9) in/lbf(IoglOof) H(9,9) in/lbf OoglO of) 





Figure 6- 1 Experimental FRF vs multiple mode matrix solutions at modes 1 and 2 using 
a 1 Hz bandwidth with a 3 point frequency samplings at each mode. 
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H(9,9) iiVlbf(loglOof) H(9,9) in/lbf (loglOof) 




Figure 6- 2 Experimental FRF vs multiple mode matrix solutions at modes 1 through 4 
using a 1 Hz bandwidth with a 3 point frequency samplings at each mode. 
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B. MULTIPLE MODE INTEGRAL SOLUTIONS 



We wish to extend single mode integral solutions to multiple modes of the system. 
As with the multiple mode matrix solutions we simply take as the path of integration a 
path which is the set theoretic union of suitable paths at each of the modes under 
consideration. As with the matrix solutions one can chose the weighting function W(s) to 
be 1/(Q) or Q to achieve lower or higher mode accuracy. As stated in the previous section 
it is required that a 'good' error set is known. Figure 6-3 is a comparison of experimental 
and uncorrected FRF versus multiple mode integral solution corrected FRF with the 
solution having been computed over modes 1 and 2 using a 1 Hz bandwidth at each mode 
with a 3 point frequency sampling over the bandwidth at each of the modes. Figure 6-4 is 
a conqjarison using a solution computed over mode 1 through 4 again using a 1 Hz 
bandwidth at each mode with a 3 point frequency sampling over the bandwidth at each of 
the modes. Figure 6-5 is a comparison plot of the multiple mode matrix and integral 
solutions over modes 1 through mode 4 using a 1 Hz bandwidth at each mode with a 3 
point frequency sampling over the bandwidth at each of the modes. 
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H(9,9) in/Ibf(loglOof) H(9,9) in/lbf (loglO of) 





Figure 6- 3 Experimental FRF vs multiple mode integral solutions at modes 1 and 2 using 
a 1 Hz bandwidth with a 3 point frequency samplings at each mode. 
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H(9,9) in/lbfaoglOof) H(9,9) in/lbf OoglOof) 





Figure 6- 4 Experimental FRF vs multiple mode integral solution at modes 1 through 4 
using a 1 Hz bandwidth with a 3 point frequency samplings at each mode. 
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H(9,9) in/Ibf(loglOof) 




Figure 6- 5 Experimental FRF vs multiple mode matrix and integral solutions at modes 1 
and 2 using a 1 Hz bandwidth with a 3 point frequency samplings at each mode. 
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C. SINGLE POINT MULTIPLE MODE SOLUTIONS 



The results of sections A and B can be performed using partitions consisting of a 
single point. Figure 6-6 shows the results of performing a multiple mode matrix solution 
over the first three modes of our spatialfy incomplete beam using a single point partition at 
each of the modes under consideration. Figure 6-7 shows the results of performing a 
multiple point integral solution over the first three modes of our spatially incomplete beam 
using a single point partition at each of the modes under consideration. Figure 6-8 
compares the multiple mode matrix and integral solution computed over the first three 
modes of our spatialfy incomplete beam with the analytic FRF. In figure 6-9 we compare a 
multiple mode matrix solution computed over the first three modes of a spatially complete 
beam using a single point partition at each of the modes with the beam's spatially complete 
analytic FRF. In figure 6-10 we compare a multiple mode integral solution computed over 
the first three modes of a spatially complete beam using a single point partition at each of 
the modes with the beam's spatially complete analytic FRF. Figure 6-1 1 con^ares multiple 
mode matrix and intgral solution computed over the first three modes of a spatialfy 
con^lete beam with the beam's spatially con^lete analytic FRF. Figures 6-9 and 6-10 
show that .the single point multiple mode solutions are exact. 
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H(9,9) in/lbfOoglOof) H(9,9) in/lbf (loglOof) 





Figure 6- 6 Experimental FRF vs multiple mode matrix solution computed at modes 1 
through 3 using a 1 Hz bandwidth with a single point frequency samplings at each mode 
for a spatially incomplete beam. 
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H(9,9) in/lbf(loglOof) H(9,9) in/lbf (loglO of) 





Figure 6- 7 Ejq)erimental FRF vs multiple mode integral solution computed at modes 1 
through 3 using a 1 Hz bandwidth with a single point frequency samplings at each mode 
for a spatially incomplete beam. 
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Figure 6- 8 Experimental FRF vs multiple mode matrix and integral solutions computed at 
modes 1 through 3 using a 1 Hz bandwidth with single point frequency samplings at each 
mode for a spatially incomplete beam. 
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H(6,6) in/lbf (logic of) H(6,6) in/lbf (loglOof) 





Figure 6- 9 Experimental FRF vs multiple mode matrix solution computed at modes 1 
through 3 using a 1 Hz bandwidth with a single point frequency samplings at each mode 
for a spatially complete beam. 
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H(6,6) in/lbf(loglOof) H(6,6) in/lbf (loglOof) 





Figure 6- 10 Experimental FRF vs multiple mode integral solution computed at modes 1 
through 3 using a 1 Hz bandwidth with a single point frequency samplings at each mode 
for a spatially complete beam. 
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H(6,6) in/lbf(loglOof) 




Figure 6-11 Experimental FRF vs multiple mode matrix and integral solution computed 
at modes 1 through 3 using a 1 Hz bandwidth with a single point frequency samplings at 
each mode for a spatially complete beam. 
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Figure 6- 12 Experimental FRF vs multiple mode matrix solutions computed over 1, 2, 3, 
and 4 modes using a 1 Hz bandwidth with single point frequency samplings at each mode 
for a spatially complete beam. 
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vn. CONCLUSIONS / RECOMMENDATIONS 



A. SUMMARY 

Frequency Domain Structural Identification using the structural synthesis 
transformation was performed on a simulated ejq>erimental fi’ee-fi'ee beam model. The 
identification was performed for both the spatially complete and the spatially incomplete 
case. 

• SST based structural identification provides an exact solution for the 
identification of FE modeling errors given spatially complete data. 

• For spatially inconplete systems SST based structural identification provides 
a fi'equency dependent solution which is not suitable for finite element modeling error 
correcting. 

• Using both matrix and integral based techniques the SST can provide single 
mode solutions which are frequency independent correction matrices AK, AM, AC 
which approximately corrects the reduced FE model within a frequency bandwidth of 
the experimental system mode under consideration. 

B. CONCLUSIONS 

This thesis has clearly demonstrated that for spatially incomplete systems, single 
mode frequency independent solutions can be found which correct the reduced finite 
element model in a neighborhood of a given mode of the ejq)erimental system. It has also 
been shown that the concept of a single mode solution can be expanded to that of multiple 
mode solutions which are frequency independent correction matrices which approximately 
corrects the reduced FE model throughout a frequency bandwidth which includes more 
than 1 mode of the experimental system. 

C. RECOMMENDATION 

The purpose of this thesis was to find frequency independent multiple mode 
solutions which could be used to approximately correct reduced finite element models. 
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Although satisfectory results were obtained, investigation is still required in the following 
areas: 

•Determine the relationship between errors in a full FE model and those of an 
associated reduced FE model. 

•Determine a method to 'pullback' reduced FE model correction to full order 
FE model corrections. 

•Investigate use of the integral formulation of reference (5) in analysis of the O- 
set system. 
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APPENDIX 



The following is a brief description of MATLAB routines employed in this thesis; 

•SST.M - Generates MAT files containing spatially complete or spatially incomplete 
experimental FRF, analytic FRF, localization matrix, and SST solutions. 

•PLOTSST.M - Uses files produced by SST.M to generate plots used in chapters 3 and 4. 

•CHAPS .M - Sets parameters and calls SST.M and PLOTSST.M to generate figures for 
chapter 3. 

•CHAP4.M - Sets parameters and calls SST.M and PLOTSST.M to generate figures for 
chapter 4. 

•CHAPS. M - Generate figures for chapter 5. 

•CHAP6.M - Generate figures for chapter 6. 

•SETUP.M - Generates MAT files for FE models. 

•BEAMMDL.M - Setup FE models. 

•FSTATIC.M - Perform static reduction of mass and stiffness matrices. 

•FERS_TAM.M - Perform IRS reduction of mass and stiflfiiess matrices. 

•FREQMODE.M - Returns FE model frequencies. 

•NDX3D.M - Indexes a series of 2-D matrices into a single 3-D matrix. 

•INTSUB.M - Decomposes impedance matrix into mass, stiffness, and damping matrices using 
integral techniques. 

•MYTRAP.M - Computes integrals using trapezoidal method 
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SST.M 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%General setup%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
clear,clc ;%%clear workspace, 

closeall ;%close aiiy open figure windows 

j=sqrt(-l) ; 

load sstcoiif ;%Ioad A-set and O-set 
aset(oset)=[]; 
aset_size=lengdi(aset); 
oset_size=Iengtli(oset); 
load beamdata 
if lengtli(oset)== 0 
complete=l; 
else 

coiiiplete=0; 

end 

%%setupplot labels%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
if (oset==[])&(posofMasseiT(l >^0)&(posofStifreiT(l 
casenanie=(’Spatially Complete Mass & StifQiess error'); 
elseif (oset=[])&(posofMasserr(l)-^0) 
casename=('Spatially Complete Mass error"); 
elseif (oset==[])&(posofStifrerr(l>^0) 
casename=('Spatially Complete Stiflness errof); 
elseif (posofMasserr( 1 )-=0)&(posofStitTerr( 1 )-=0) 
casename=(’Spatially Incomplete Mass & Stiffness error'); 
elseif (posofMasserrf 1 )^0) 
casename=('Spatially Incomplete Mass error*); 
elseif (posofStifferr( 1 )-^0) 
casename=('Spatially Incomplete Stiffness error*); 
end 

%%get reduced A set O set frequencies%%%%%%%%%%%%% 

kexto=k_anal(oset,oset); %% stiffness 
mexto=m__anal(oset,oset); %% and mass matrices 
if complete 
kstat=k_anal; 
mstat=m_anal; 
else 

if static == 0 

[kstat,nistat]=fstatic(k_anal,m_aiial,oset,aset);%% get reduced K & M 
elseif static == 1 

[kstat,mstat]=firs_lam(k_anal,m_anal,oset,aset);%% get reduced K & M 

else 

kstat=k_anal(aset,aset); 

instat=m_anal(aset,aset); 

end 

end 

cstat=sqrl(-l )*struc_damping. ^kstat; 

[u,lambdaa,c]=freqinode(k_anal,m_anal) ;%get freqs 
[u,lambdared,c]=freqmode(kstat 4 nstat) ;% 

[u,lambdax,c]=ffeqmode(k__exp,m_exp) ;% 

[u,lambdaexto,c]=freqmode(kexto,mexto);% 

omegaa=sqrt( lambdaa ); 

omegax=sqrt(lainbdax); 

omega red= sqrt( lambdared); 

omegaexto=sqrt( lambdaexto ); 

for i=l .4 % diis to delete fixed bod}’ modes 
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SST.M 



if real(omegax( 1 )) < 1 0^(-3 ) 

omegax=omegax(2:lengUi(oniegax)); 

end 

if real(omegaa( 1 )) < 
omegaa=oniegaa(2:leiigtli(omegaa)); 
end 

if complete==0 

if real(omegared( 1 )) < 10^(-3) 
omegared=omegared(2 : lengtli(oniegared)); 
end 

if real(omegaexto( 1 )) < 1 0^(-3) 
omegaexlo=omegaexto(2:lengUi(omegaexto)); 
end 

end 

end 

%%%%%tmd true mass and stifliiess errors%%%%%%%%%%%%%%%%%%% 
true_stiflhess=k_exp-k_aual; 
stifliiess_cset=fiiid(diag(true__stifliiess)); 
true_mass=m_exp-m_aiial; 
mass_cset=find(diag( true_mciss )); 
true_dainping=c_exp-c_anal; 
damping_cset=find(diag(true_damping)); 
if lengtli(mass_cset) > 0 
for i=l:length(inass_cset) 

x=find(stiflhess_cset==niass_cset(i)); 
if~(lengtli(x)> 0) 

stiffiiess_cset= [ stiffness_cset; inass_cset(i )] ; 
end 
end 
end 

if lenglli(damping_cset) > 0 
for i=l.ieiigtli(daniping_cset) 

x=find(stiflhess_cset==dainping_cset(i)); 
if -<lengdi(x) > 0) 

stiffiiess_cset=[sliffiiess_cset; damping__cset(i)]; 
end 

end 

end 

tnie_cset=sort( stiffiiess_cset) 

save true_errs true_cset true_mass true_stiffiiess true_daniping 
clear true_mass true_stifiiiess tnie_damping 

clear area eeii force g gc gcx gk gkx gin guLX goblc goblcx goblk goblkx 
clear goblin goblnix ke lainbdaa lanibdaexto lambdared lainbdax lumpdainp 
clear lunipmass lumpspring me pho posofDaniperr posofMasserr posolStifferr 
clear stiffiiess_cset 

%%%setup frequency range for aiialysis%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

freqtop=oniegax(liiglimode)+. 1 *oinegax(liiglunode); 

freqbottom=omegax(lowmode)-. 1 *omegax( lowmode); 

freqbottoiu= 10*2 *pi; %start at 1 0 liz 

freqtop=1200*2*pi; %end at 1200 liz 

%%number of points to plot 

numpoints=400;%%%%%SET NUMBER OF POINTS TO USE FOR THE SIMULATION 
fmeness=numpoints; 

vv=freqbottom:(freqtop-freqbottom)/(numpoints- 1 ):freqtop; 

%%Ixcate C set%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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SST.M 



keepgomg=’n' 

\vfi-eq=omegax(lo\vinode) ; 
nextiiiode=lowniode; 
tol= 1 ; 

defaultstepsize=. 1 ; 
while (keepgoiiig=’n’) 

z_anal_red=kstat+j ♦wfreq*cslat-wfreq^2 *mstat ; 

h_anal_red=iiiv( z_anal_red) 

z_exp=k_exp+j * w£req*c_exp-wlreq^2 *in_exp ; 

h_exp=inv(z_exp) ; 

h_exp=h_exp(aset,aset) 

L=z_aiial_red*(h_anal_red-h_exp) *z_anal_red ; 

ci=fmd(abs(diag(L))>tol) ; 
lemp_cset_size=lengtli(ci); 
if length(ci) >= 1 

cset_size=lengtli(ci) ; 

figure(l) 

plot(aset, abs(diag(L))) ; 
hold ou 

plot(aset(ci),zeros( 1 ,lengtli(ci)),’x') 
tallest=max(abs(diag(L))) ; 

title([TOL = ’4ium2str(tol),' Cset size = ’,mt2slr(cset_si2e),’ Omega = '^ium2str(wtreq/(2*pi)),' Hz’]) 
ylabel('lbf7in’) 
if complete 
xlabel(’DOF) 
else 

xlabel(’ASETDOF) 

end 

hold off 

k=iueuu( ' Choose an action 
'hicrease Tolerance . 

Decrease Tolerance 
Increase frequency 
Decrease frequency 
' Previous mode 
' Next Mode 
'Set Print switch on',... 

' hicrease stepsize ',... 

' Decrease stepsize 
• Go ’); 

ifk==l 

stepsize=defaultstepsize; 

ctr=0; 

wliile lengtli(ci) >= teinp_cset_size 
ctr=ctrfl; 
if ctr==10 

stepsize=stepsize* 1 0, 
cti=0 ; 
end 

tol=min(abs(tol+stepsize^(tol)),tallest) ; 
c i=fmd(abs(diag(L))>tol ) 
iflength(ci)>= 1 
cset_size=lenglh(ci) ; 

figure(l) 

plot(aset, abs(diag(L))) ; 
hold on 

plot(aset(ci),zeros( 1 ,lengtli(ci)),'x’) 
tallest=niax(abs(diag(L))) ; 

title([TOL = ',nuni2slr(lol),' Cset size = ',iiit2str(cset_size),’ Omega = ',nuiu2str(\vfreq/(2*pi)).' Hz']) 
ylabel('lbfyin') 
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SST.M 



if complete 

xlabel(’DOF') 

else 

xlabelCASET DOF) 

end 

hold off 
end 
end 

elseif k==2 

stepsize=defaullstepsize; 

ctr=0; 

while length(ci) <= tenip_cset_size 
ctr=ctrH; 
if cti==10; 
stepsize=stepsize*2; 
ctr=0; 
end 

tol=max(abs(tol-stepsize*(tol)),.0000 1 ) ; 
ci=find(abs(diag(L))>tol) 
if lengtli(ci)>= 1 

cset_si 2 e=leiigtli(ci) ; 

figure(l) 

plot(aset, abs(diag(L))) ; 
hold on 

plot(aset(ci),zeros( 1 ,leiigtli(ci)),'x’) 
lallest=max(abs(diag(L))) ; 

title([TOL = ’,num2str(tol),' Cset size = *,int 2 str(cset_si 2 e),’ Omega = ',iium2str(v\'treq/(2*pi)),' 

ylabel('lbf7in’) 
if complete 
xlabel(DOF') 
else 

xlabelCASET DOF) 
end 

hold off 
end 
end 

elseif k=3 

while lengtli(ci) == temp_cset_size 

wfreq=mm( wfreq+.OOOO 1 *(Creqtop-wfreq),freqtop); 

z_anal_red=kstat+j *w(count)*cstat-wlreq'^2*mstat ; 

h_anal_red=inv(z_anal_red) 

z_exp=k_exp+j * wfreq*c_exp-\vfreq^2 *m_exp ; 

h_exp=inv(z_exp) ; 

h_exp=h_exp(aset,aset) 

L=z_anal__red *(h_anal_red-h_exp ) *z_anal__red 
tallest=max(abs(diag(L))) 
tol=inin(tol+.001 *tallest,tallest) ; 
ci=fmd(abs(diag(L))>tol) ; 
cset_size=length(ci) ; 

figiire(l) 

plot(aset, abs(diag(L))) ; 
hold on 

plot(aset(ci),zeros( 1 ,lengtli(ci)),V) 

title([TOL = ',num2slr(tol),’ Cset size = ',int2slr(cset_size),' Omega = ',nuni2str(wfreq/(2*pi)),’ Hz’]) 

ylabel(’lbf7in’) 

if complete 

xlabel(’DOF) 

else 

xlabel(ASET DOF) 
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end 

hold off 
end 

elseif k==4 

wfreq=max(vv1req-.05*(wfreq-freqbottom),freqbottoin); 
elseif k==5 

nextmode=max(lowraode,nextmode-l ); 
w£req=oinegax(nextinode); 
elseif k==6 

nextmode=min(liiglimode,nextmode+ 1 ); 

\vfreq=omegax( nextmode); 
elseif k==7 

pswitcli=y ; 

elseif k==8 

defaultstepsi 2 e=defaultstepsi 2 e/ 10 ; 
elseif k==9 

defaultstepsize=defaullstepsize* 1 0; 

else 

keepgoing=y 
cset=aset(ci); 
cset_rel=ci ; 

end 
else 

tol=tol*.9; 

end 

end 

close(l) 

save extrtset cset cset_rel 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%% 

%%%%%%%%%%%%%%%%%%FORCETHECSET%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%csel=[79 11 13] %%%%%%%%THE INCOMPLETE CASE CSET%%%%%%%%%%%% 

%cseLrel=[4 5 6 7]%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%«/c%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%% 

%csel=[4 5 6 7 8 9 10 1 1]%%%%%%%%%%%%%THE COMPLETE CASE CSET%%% 

%cset_.rel=[4 5 6 7 8 9 10 11]%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%% 

cset_size= length( cset); 
save L L 
clear L 

clear kexta kexto mexta iiiexlo 
pack 

skyfull=zeros(nimidof,ntundoO; 
coimt=l ; 

for index 1 = 1 :numdof 
for index2=indexl:numdof 

sky full( index 1 ,index2 )=count; 
skyfull(index2,indexl )=count; 
count=count+l; 
end 
end 

skyred=zeros(aset_size,aset_si 2 e); 
count= 1 ; 

for index 1 = 1 :aset_size 
for index2=indexl :aset_size 

skyred( index 1 ,index2 )=count; 
skyred( mdex2 index 1 )=count. 
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couiit=couiit+l; 

end 

end 

skycset=zeros(cset_size,cset_size); 
count=l ; 

for index 1 = 1 :cset_size 
for index2=indexl:cset_size 

skycset(index 1 ,index2)=couiil; 
skycset(index2, index 1 )=couiit; 
coimt=comit+l; 
end 
end 

full_hoideT=zeros( 1 ,(mmidoP(nmndof+ 1 ))/2); 
red_holder=zeros( 1 ,(aset_size*(aset_size+* 1 ))/2 ), 
cset_holdei=zeros( 1 ,(cset_size’(cset_size^l ))/2); 

if meters 

waitbar_handle=waitbar(0,'Computing Experimental FRF'); 
else 

dispCGetting experimental FRF') 
end 

H_EXP=Q; 

for coimt=l .nmnpoints 
if meters 

waitbar(count/nunipoints); 

end 

z_exp=k_exp+j * w(comit)*c_exp-w(comU)^2 *m_exp; 

h_exp=inv(z_exp); 

h_exp=h_exp(aset,aset); 

skyindex=l; 

for index 1 = 1 :aset_size 

for index2=index 1 :aset_size 
red__Iiolder(skyindex)=h_exp(indexl,index2), 
skyindex=skyindex+ 1 ; 
end 
end 

H_EXP=[H_EXP redjiolder’]; 
end 

if meters 

cIose(waitbarJiandle) 

end 

saveH_EXPH_EXP 
clear H_EXP 
pack 

if meters 

waitbarJiandle=waitbar(0,'Computing Experimental Impedence'); 
else 

dispCGetting experimental Impedence') 
end 

Z_EXP=[]; 

for count= 1 :nmnpoints 
if meters 

waitbar(count/numpoints); 

end 

z_exp=k_exp+j ♦\v(count)*c_exp-\v(coiint)'^2 *m_exp; 

h_exp=inv(z_exp); 

h_exp=h_exp(aset,aset); 

z_e\p=inv(h_exp); 
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skyindex=l ; 

for indcxl=l:asel_size 

for iiidex2=indexl;aset_size 
red_holder(sk 7 iiidex)=z_exp(index 1 ,index2); 
skyiiidex=skyindex+ 1 ; 
end 
end 

Z^EXP=[Z_EXP redjiolder*]; 
end 

if meters 

close(\vailbar_handle) 

end 

saveZ^EXPZ^EXP 
clear Z_EXP 
pack 

if meters 

waitbar_han(Ue=wailbar(0,’Conipuling Reduced Analytic FRF'); 
else 

disp(’Getting Reduced FRF’) 
end 

H_ANAL_RED=[]; 
for count=l:iiumpoints 
if meters 

wailbaifcount/numpoinls); 

end 

z_anal_red=kstat+j * w(count)*cslat-w(count)'^2 *mstat; 

h_anal_red=inv(z_anal_red); 

skyindex=l ; 

for index 1 = 1 :asel_size 

for index2=indexl :aset_size 
red__holder(skyindex)=h_anal_red( index 1 ,index2); 
skyindex=skyiiidex+l ; 
end 
end 

H_ANAL^RED=[H^ANAL_RED redjiolder*]; 
end 

if meters 

close(wailbarjiandle) 

end 

save H^ANAL_R H_ANAL_RED 

clear Z_ANAL_RED H_ANAL_RED kexta kexto mexta mexto 
pack 

%% compute and save L matrix diagonal as a function of frequency%%%%%%%% 
if meters 

waitbar Jiandle=wailbar(0,'Compuling L Diagonals'); 
else 

disp('Gelting L Matrix') 
end 

L^DIAGS=[]; 
for couiU= 1 inmiipoints 
if meters 

\vaitbar(count/numpoints); 

end 

z_anal=k_anal+j*w(counl)*c_anal-w(coiml)'^2*m_anal; 
z_anal_red=kstat+j * w(count)*cslal-\v'(coiuil)^2 *mslat; 
h_anal_red=inv(z_anal_red); 
z_exp=k_expH'j *c_exp-\v (count)^2 *in_exp; 
h_exp-inv(z_exp); 
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h__exp=h_exp(aset,aset); 

L=z_anal_red*(h_anal_red-h_exp)*z_aiial_red; 

temp_l__diags=diag(L); 

L_DIAGS=[L_DIAGS temp_l_diags(;)]; 
end 

if meters 

close(waitbar_handle) 

end 

save L^DIAGS L_DIAGS 
clear L_DIAGS z__aiial temp_l_diags L 
clear z_exp z_anal_red u v h_exp h_anal_red 
pack 

if meters 

waitbar_handle=waitbar(0,'Computmg Z Analytic'); 
else 

disp('Gettiiig Z anal') 
end 

Z^ANAL=0; 

for coimt= 1 :numpoints 
if meters 

waitbar(count/numpoints); 

end 

z__anal=k_anal+j * w(count)*c_anal-w(count)^2 *m_anal; 

skyindex=l ; 

for indexl=l:nunidof 

for index2=indexl;numdof 
full_holder(skyindex)=z_anal(index 1 ,index2 ); 
skyindex=skyindex+ 1 ; 
end 

end 

Z_ANAL=[Z_ANAL fulLholder*]; 
end 

if meters 

close(waitbarJiandle) 

end 

save Z_ANAL Z_ANAL 
clear Z_ANAL z_anal 
pack 
if meters 

\vaitbar_handle=waitbar(0,'Computing Reduced Z Anal}!^'); 
else 

dispCGetting Reduced Z anal') 
end 

Z_ANAL_RED=[]; 
for count= 1 .numpoints 
if meters 

waitbar(count/numpomts); 

end 

z_anal_red=kstat+j * w(count) * cstat-\v(count)^2 * mstat; 

skyindex=l ; 

for indexl = 1 :aset_size 

for index2=indexl :aset_size 
red_holder(skyindex)=z_anal_red( index 1 ,index2 ); 
skyindex=skyindex+ 1 ; 
end 

end 

Z_ANAL_RED=[Z^ANAL_RED redjiolder^]; 
end 

if meters 
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close(waitbarjiandle) 

end 

save Z_ANAL_R Z_ANAL_RED 
clear Z_ANAL_RED z_anal_red 
pack 

%%%compute DZ%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
if meters 

waitbar_handle=waitbar(0,'Computing DZ. . 
else 

disp(’Getting DZ’) 
end 
DZ=[]; 

for count=I:numpoints 
if meters 

waitbar(count/nuinpoints); 

end 

z_anal_red=kstat+j ♦ w(count)’^cstat-\v(count)"^2 *mstat; 

h_anal_red=inv(z_anal_red); 

z_exp=k_exp+j*w(count)*c_exp-\v(count)^2*m_exp; 

h__exp=inv(z_exp); 

h_exp=h__exp(aset,aset); 

hacc=h_anal_red(cset_rel,cset_rel); 

lixcc=h_exp(cset_rel ,cset_rel); 

dz=inv(inv(inv(hacc)*(liacc-hxcc)*inv(hacc))-hacc); 

skyindex=l ; 

for index 1 = I :cset_size 

for index2=index 1 :cset_size 
csetholdeif skyindex)=dz(mdex 1 ,index2 ); 
skyindex=skyindex+ 1 ; 
end 

end 

DZ=[DZ csetjiolder*]; 
end 

if meters 

close(\vaitbar_handle ) 
end 

save DZ DZ 

clear hacc lixcc dz h_exp z_anal_red h_anal_red z_exp 
pack 

%%%decompose DZ into DM, DK, & DC%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
if meters 

waitbarJiandle=\vaitbar(0,T)econiposing DZ...’); 
else 

disp(’Decomposmg DZ’) 
end 
DK=(]; 

DM=Q; 

DC=n; 

dzl =zeros(cset_size,cset_size); 
dz2=zeros(cset_size,cset_size); 
dz3=zeros(cset_si2e,cset_size); 
for i=l :numpoints-2 
if meters 

waitbar(i/(numpoints-2)) 

end 

dztemp=DZ(ndx3d([l cset_size*(cset_size+ 1 )/2 nmnpomts], 1 , 1 :cset_size^(cset_size+ 1 )/2,i)); 

skyiiidex=l; 

for uidex 1 = 1 :csct size 
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for index2 = index 1 :cset_size 
dz 1 ( index 1 ,index2)=dztemp( skyindex ); 
dzl (index2,iiidex 1 )=dztemp(skyindex); 
skyindex=skyindex+ 1 ; 
end 
end 

dztemp=DZ(ndx3d([ 1 cset__size*(cset__size+l )/2 luunpoints], 1 , 1 :cset__size*(cset_size^- 1 )/2,i+l )); 

sk>index=l; 

for index 1 = 1 :cset_size 

for index2=indexl:cset_size 
dz2(index 1 ,index2 )=dztemp(skyindex); 
dz2(index2,index 1 )=dztemp( skyindex); 
skyindex=skyindex+ 1 ; 
end 
end 

dztemp=DZ(ndx3d([l cset_size*(cset_size+ 1 )/2 numpoints], 1 , 1 :cset_sizc*(cset_size+ 1 )/2,i+2)); 

skyindex=l; 

for index 1 = 1 :cset_size 

for index2=indexl:cset_size 
dz3(index 1 ,index2)=dztemp( skyindex); 
dz3( index2 ,index 1 )=dztemp( skyindex ); 
skyindex=skyindex+ 1 ; 
end 
end 

temp=[eye(cset_size) -w(i)^2*eye(cset_size) +j*w(i)*eye(cset_size); 
eye(cset__size) -w(i+l)^2*eye(cset_size) +j’^w(i+l)*eye(cset_size); 
eye(cset_size) -’w(i-i-2)^2*eye(cset_size) +j *vv(i+2)*eye(cset_size)]\[dzl;dz2;dz3]; 
ktemp=temp( 1 :cset_size,:); 
mtemp=temp(cset_size+ 1 :2 *cset_3ize,: ); 
ctemp=temp(2 * cset__size+ 1 :3 * cset_size,: ); 

skyindex- 1 ; 

for index 1 = 1 :cset_size 

for index2=indexl:cset_size 
cset_holder(skyindex)=ktemp(indexl,index2); 
skyindex=skyindex+ 1 ; 
end 
end 

DK=pKcsetJiolder(:)]; 

skyindex=l 

for indexl = 1 :cset_size 

for index2= index l:cset_size 
cset Jiolder(skyindex)=mtemp(index 1 ,i ndex2 ); 
skyindex=skyindex+l ; 
end 
end 

DM=(DM cset_holder(:)]; 

skyindex=l ; 
for index 1 = 1 :cset_size 

for index2=indexl :cset_size 
cseMiolder(skyindex)=ctemp(index 1 ,index2); 
skyindex=skyindex+ 1 ; 
end 
end 

DC=[DC cset_holder(:)]; 
end 

if meters 

close(vvaitbar_handle) 
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end 

save DM DM 

save DK DK 

save DC DC 

clear DC DK DM 

save exp c_exp k_exp ra_exp 

clear c_exp k__exp m_exp 

save stat kstat mstat cstat 

clear kstat mstat cstat 

save anal k_aiial m_anal c_anal 

clear k_aiial m_anal c_anal 

pack 

if complete 
load Z__ANAL 
Z=Z_ANAL; 
clear Z_ANAL 
else 

load Z_ANAL_R 
Z=Z_ANAL_RED; 
clear Z_ANAL_RED 
end 

CORRHA=[]; 

tempza=zeros(aset_si 2 e,aset_size); 
teinpz=zeros(cset_size,cset__size); 
if meters 

waitbar_handle=waitbar(0,’lnstalling DZ. . 
else 

disp(lnstallmg DZ') 
end 

for i”l:numpoints 
if meters 

waitbar(i/nunipoints) 

end 

ztemp=Z(ndx3d([l aset_size*(aset_size+l)/2 (numpoiiits-(i-l ))],1,1 :a$et_size*(aset_size+l)/2,l)); 

skyindex=l; 

for index 1 = 1 :aset_size 

for index2=indexl:aset_size 
tempza(indexl,index2)=ztemp( sky index); 
tempza(index2, index 1 )=ztemp( sky index); 
skyindex=skyindex+ 1 ; 
end 
end 

Z(:,iHl; 

ztemp=DZ(ndx3d([ 1 cset_size*(cset__size+l )/2 (numpoints-(i-l ))], 1 , 1 :cset_size*(cset_size^- 1 )/2, 1 )); 

skyindex=l; 

for index 1 = 1 :cset_size 

for index2=indexl:cset_size 
tempz(index 1 ,iiidex2)=ztemp(skyindex); 
tempz(index2,index 1 )=ztemp(skyindex); 
skyindex=skyindex+ 1 ; 
end 
end 

tenipza(cset_rel,cset_rel)=tempza(cset_rel,cset_rel>i-tempz; 

DZ(:,1)=[], 

tempha=inv(tempza); 

CORRHA=[CORJRHA templia(:)]; 
end 

if meters 

close(waitbar_handle) 
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end ' 
clear DZ Z 

save CORRHA CORRHA 
clear CORRHA 

load DK; 
load DM; 
load DC; 

load slat 

CORRHAD=0; 

tempza=zeros(aset_size,aset_size); 

tempz=zeros(cset_size,cset_size); 

tempk=zeros(cset__size,cset__size); 

tempm=zeros(cset__size,cset__size); 

tempc=zeros(cset_size,cset_size); 

if meters 

waitbar_handle=\vaitbar(0,'lnstalliiig DK, DM, & DC...'); 
else 

dispCInstalling DK, DM, DC) 
end 

for i= 1 :iuunpoiiits-2 
if meters 

waitbar(i/(numpoiiits-2)) 

end 

kcorrected=kstat; 

mcoirected=mstat; 

ccorrected=cstat; 

temp=DK(ndx3d([l cset_size*(cset_size+l )/2 (numpoints-(2))], 1 , 1 :cset__size*(cset__size+ 1 )/2,l )); 

skyindex=l; 

for index 1 = 1 :cset_size 

for index2=indexl:cset_size 
tempk(indexl,index2)=temp(skyiiidex); 
tempk(index2, index 1 )=temp(skyiudex); 
skyindex=skyindex+ 1 ; 
end 
end 

DK(:,1)=[]; 

temp=DM(ndx3d([ 1 cset_size^(cset_size+ 1 )/2 (uimipoints-(2))], 1 , 1 :cset__size*(cset_size+ 1 )/2, 1 )); 

skyindex=l; 

for index 1 = 1 :cset_size 

for index2=indexl :cset_size 
tenipm(indexl,index 2 )=temp(sk 7 index); 
tempin( index2,index 1 )=temp(skyindex); 
skyindex=skyindex+ 1 ; 
end 
end 

DM(:,1H]; 

temp=DC(ndx3d([ 1 cset_size*(cset_size+ 1 )/2 (numpoints-(2))], 1 , 1 :cset_size*(cset_size+ 1 )/2, 1 )); 

skyindex=l; 

for index 1=1 :cset_size 

for index2=uidexl:cset_size 
tempc(index 1 ,index2)=temp(sky index); 
tempc(index2,index 1 )=temp(skyindex); 
skyindex=skyindex+l ; 
end 
end 
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DC(:,1H]; 

kcoirected(cset_rel,cset_rel)=kcoiTected(cset__rel,cset_rel)^'tempk; 

mcoirected(cset_rel,cset_rel)=incoiTccted(cset_rel,cset_rel)+teinpiii; 

ccorrected(cset_rel,cset_rel)=ccorrected(cset_rel,cset_rel)+teinpc; 

tempza=kcorrected+j * w(i+ 1 )*ccoirected-w( i+ 1 )^2 "^mcorrectcd; 
tempha=inv(tempza); 

COmiHAI>[CORRHAD templia(:)]; 
end 

if meters 

close(vvaitbar_handle) 

end 

clear teinpk tempin tempo tempza tempha 
clear kcorrected mcorrected ccorrected 
save CORRHAD CORRHAD 
clear CORRHAD DK DM DC 

%%%save tlie 

Workspace%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
save INT 
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clear 
closeall 
load INT 
titles=0; 
if psvvitch==y 
wliitebgCwliite') 
close 
end 
h=0; 

if complete 
fignum=l; 
else 

figiumi=21; 

end 

%% \vrite frequencies and system diag to diary file prt.out%%%%%%%%%%%%%%%%% 
if exist(TigOOO.out') 
delete GgOOO.out 
end 

diary figOOO.out 
fprintf(’ \n') 
fprintfC \n’) 
fpriutf('A set\iV) 
for index=l:lengtli(aset) 

fprintf([blanks(4-lengtli(int2str(aset(index)))),int2str(aset( index))]) 
if rem( index, 1 2)==0 
fprintfCV) 
end 
end 

fprintf(V) 
fprintfC\n') 
fprintf('0 setW) 
for index=l:lengtli(oset) 

fprintf([blanks(4-lengtli(int2str(oset(index)))),int2str(oset(index))]) 
if rem( index, 12)==0 
fprintflV) 
end 
end 

fprintf(V) 

fprintfC\n') 

fprintf(’Computed C setW) 
for index=l:lengtli(cset) 

fprintf([blanks(4-lengtli(int2str(cset(index)))),int2slr(cset(index))]) 
if rem( index, 1 2)==0 
fprintf(V) 
end 
end 

fprintfC\n') 

fprintffV) 

fprintf(True C setW) 

for index=l:lengtlr(true_cset) 

fprintf([blanks(4-lengtli(int2str(true_cset(index)))),int2str(true_cset(index))]) 
if rem( index, 12)==0 
fprintf(V) 
end 
end 

fprintf(^^l’) 

fprintf(W) 

bigger=inax(lengtli(oinegaa),lengtli(omegax)); 
if complete 
allfreqs=[ ... 
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[omegaa/(2*pi); zeros(bigger-lengtli(oniegaa),l)] ... 

[omegax/(2*pi); zeros(bigger-leiigtli(omegax), 1 )] ... 

]; 

else 

allfreqs=[ ... 

[oniegaa/(2*pi); zeros(bigger-lengtli(omegaa), 1 ) ] . . . 

[omegax/(2*pi); zeros(bigger-length(omegax),l) ]... 

[omegared/(2*pi); zeros(bigger-length(omegared),l) ] ... 

[omegaexto/(2*pi); zeros(bigger-lengtli(oinegaexto), 1 )]. . . 

]; 

end 

[aaa bbb]=size(all{reqs); 
fprintf(’System Frequencies (Hz)\n') 
if complete 

fprintf(’ Anal Exp\n’) 
bbb=2; 
else 

fprintfC Anal Exp Red Oset\n') 
bbb=4; 
end 

for indexl=l:aaa 
prlline=[]; 
for index2=l:bbb 
if allfreqs(indexl ,index2)==0 
prtfreq=’ 
else 

prtfreq=(blaiiks(8»lengtli(sprintf(%4g',allfreqs(index 1 ,iiidex2 )))),. . . 
sprintfC%4g*,allfreqs( index 1 ,index2))]; 
end 

prtline=[prtliiie, prtlreqj; 
end 

prtliiie=[prtline, ^n^; 
fprintf(prtline) 
end 

diary off 

%dos(*type figOOO.out > Iptl &’); 

xx=omegax(fmd(omegaa(lovvinode) <= oniegax & omegax <= omegaa(liiglunode)))/(2*pi); 

>-x 1 =ones(length(xx), 1 ); 

aa=omegaa(fmd(omegaa(lowinode) <= omegaa & omegaa <= oniegaa(higlunode)))/(2*pi); 
ya 1 =ones(lengtli(aa), 1 ); 

red=omegared(find(oniegaa(lowmode) <= omegared & omegared <= oniegaa(higlunode)))/(2*pi); 
yred 1 =ones(leiigtli(red), 1 ); 

e\1o=omegaexto(fnid(omegaa(lo\vmode) <= omegaexlo & omegaexlo <= omegaa(liighniode)))/(2*pi); 
yexto 1 =ones(lengtli(exto), 1 ); 
if complete 

x0=max([xx(lengtli(xx)) aa(lengtli(aa))] ); 
else 

x0=max([xx(length(xx)) aa(lengtli(aa)) red(length(red)) exto(lengtli(exto))]); 
end 

clear oniegax omegaa omegared omegaexto 

clear L Z_ANAL c_anal c__exp conn h_exp h_anal_red k_anal k_exp 
clear kexta kstat kexto m_anal m__exp mexto instat temp_l_diags 
if complete 

mid_index=round(lengtli(cset)/2); 

else 

mid_index=round(lengtli(cset)/2); 

end 

load H_ANAL_R 
load H_EXP 
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h=figure(h+l); 
figiiuni=fignuiii+l ; 

plol(w/(2*pi),logl0(abs(H_ANAL_R£D(niix3d([l aset_size*(aset_size+l )/2 numix)ints],... 
I ,sk>Ted(cset_reI(niid_index),cset_rel(mid_iudex)),V)))),'r-’,. . . 
\v/(2*pi),loglO(abs(H_EXP(iidx3d([l aset_size*(aset_size+ 1 )/2 iuunix>ints],... 
l,skyred(cset_rel(iuid_iiidex),cset_rel(niidjndex)),’:')))),’b-.’) 
v=axis; 

ylabel(['H(',mt2str(cset(niid_mdex)),V,mt2slr(cset(niid_iiidex)),') iii/lbf (LoglO oO’]) 
xiabel('Oniega (Hz)') 
if titles 

title('Aiialytic FRF vs Experimental FRF') 
end 

hold on 

if abs(v(4)-v(3)) <= 1 
v2(l)=v(l); 
v2(2)=v(2); 

v2(3)=v(3)-.l*abs(v(4)); 
v2(4)=v(4)+. 1 *abs(v(4)); 
axis(v2); 
end 
grid on 

Wi=slegend('ANALYTIC','EXPERIMENTAL'); 
axes(lili); 
hold off 
if pswitch=='y' 
prtfig(fignuin) 
delete(h) 
end 

clear H_ANAL_RED 
clear H_EXP 

if complete 
h=figure(h+l); 
figiiuni=fignum+ 1 ; 
else 

h=figure(h+2); 
fignum=fignum+2; 
end 
load L 

plot(aset,abs(diag(L)),aset,abs(diag(L)),'*') 
ylabelC L(DOF) Ibf/in') 
if complete 
xlabel(’DOF') 
else 

xlabel(’ASETDOF') 
end 
grid on 
if titles 

title(['Localization Matrix Diagonal at Omega = ',iuun2str(\vfreq/(2*pi)), ' Idz']); 
end 

if pswitch==’y' 
prtfig(fignum) 
delete(h) 
end 
clear L 

h=figure(h+l); 
fignum=fignum+ 1 ; 
load L_DIAGS 
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mesh(w/(2*pi),aset,loglO(abs(L_DIAGS))) 
if titles 

titIe('Frequeiicy Dependence of Localization Matrix Diagonals') 
end 

xlabel('Oiiiega (Hz)’) 
ylabelCDOF) 

zlabelCL(DOF, omega) IbD^in (Log 10 of)') 
grid on 

if pswitch=y 
prtfig(fignum) 
delete(h) 
end 

h=figure(h+l); 
figiiuni=£igmim+l ; 
subplot(2,l,l) 

plot(w/(2*pi)Jog 10(L_DIAGS(ndx3d([aset_size 1 numpoints],cset_rel(mid_index+2), 1 ,’:’)))) 
ylabel(['Error Coord L(',int2str(cset(mid_index+2)),',',int2str(cset(inid_index+2)),’)']) 
if titles 

title('lbf7in (Log 10 of)' ); 
end 
grid on 
vl=axis; 
subplot(2,l,2) 

if lengtii(cset_rel) < length(aset) 
holdset= 1 ;lengtli(aset) ; 
holdset(cset_rel)=[] ; 
index=round(lengtli(holdset)/2); 

plot(w/(2*pi)JoglO(L_DlAGS(ndx3d([aset_size 1 numpoints],holdset( index),!, ':')))) 
ylabel(['Non-Error Coord L(',int2str(aset(holdset(index))),',’,int2str(aset(holdset( index))),')']) 
if titles 

title(Trequency Dependence of L Matrix Non-Error set Diagonal Elements' ) 
end 
v=axis; 
hold on 
v2(l)=v(l); 
v2(2)=v(2); 

v2(3)=v(3)+((v(4)-v(3))*.5Hvl(4)-vl(3))*.5; 
v2(4)=v(3>K(v(4)-v(3))*.5)+(vl(4)-vl(3))*.5; 
axis(v2); 
grid on 
end 

xlabel(’Omega (Hz)’) 
hold off 
if pswitch=='y' 
prtfig(fignimi) 
delete(h) 
end 

clear L_D1AGS 

load Z^ANAL_R 
loadZ^EXP 
h=figure(h+l); 
fignuin=fignuni+ 1 ; 

pIot(w/(2*pi),log 10(Z_ANAL_RED(ndx3d([ 1 aset_size*(aset_size+ 1 )/2 
luuii points], 1 ,skyred(cset_rel(mid_index),cset_rel(inid_iiidex)),':'))),’r~',. . . 

w/(2*pi),logl0(Z_EXP(ndx3d([l aset_size^(aset_size-f-l )/2 
iuuiipoints],l,sk>Ted(cset_reI(niid Judex), cset_rel(mid_index)), ':’))), 'b-.') 
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xlabel(’Omega'^2 (Hz)') 

ylabel(['Z(',iiit2str(cset(inid_mdex)),V,iiit2str(cset(iuid_iiidex)),') lbf7iii (LoglO oO']) 
if titles 

title('Analytic Impedance vs Experimental Impedance') 
end 
v=axis; 
hold on 

if abs(v(4)-v(3)) <= 1 
v2(l)=v(l); 
v2(2)=v(2); 

v2(3)=v(3)-. 1 *abs(v(4)); 
v2(4)=v(4)f . 1 *abs(v(4)); 
axis(v2); 
end 
grid on 

ifcoinplete== 0 

lili=slegendCANALYTIC IMPEDANCE','EXPERIMENTAL IMPEDANCE'); 
else 

Wi=slegend('ANALYnC IMPEDANCE','EXPERIMENTAL IMPEDANCE'); 
end 

axes(lili) ; 

hold off 
if p$witch==’y' 
prtfig(fignuin) 
delete(h) 
end 

clear Z_ANAL^RED Z_EXP 

load true_errs 
load DK 
h=figure(h+l) ; 
fignuin=fignum+l ; 
if complete 

plot(w(l:nmnpoints-2)/(2*pi),DK(ndx3d([l cset_size*(cset_size+l)/2 numpomts-2],l,skycset(midJiidex,mid_index),':')),'r- 

f 

w( 1 :numpoints-2)/(2*pi),true_stiniiess(cset(midJndex),cset(mid_index)). *ones( 1 ,numpoints-2),’b-.’) 
ylabel([DK(',int2stitcset(mid_index)),',’,int2str(cset(inid_index)),') lbI7in’]) 
xlabel('Omega (Hz)') 
else 

plot(\\'( 1 :numpoints-2 )/(2 *pi),log 1 0(abs(DK(ndx3d([ 1 cset_size*(cset_size+ 1 )/2 nuinpoints- 
2],1 ,skycset(inid__index,mid_index),':')))),'r~',. .. 

w( 1 :nuinpoints-2)/(2*pi),logl0((lrue_stiffness(cset(mid_iiidex),cset(inid_index)) == 0)+,.. 
true_stiffness(cset(mid_index),cset(niid_index))). *ones( 1 ,munpoints-2 ),'b-.') 
ylabel(['DK(',int2str(cset(inid_mdex)),',’,mt2str(cset(inid_mdex)),') lbf7in (LoglO ol7]) 
xlabel('Oinega (Hz)') 
end 

v=axis ; 
if titles 

title('Computed Stifiiiess Error vs True Stiflhess Error*) 
end 

v2(l)=v(l); 
v2(2)=v(2); 
if complete 
if abs(v(4)-v(3)) < 1 
v2(3)=v(3)-100*abs(v(4).v(3)); 

v2(4>=v(4>fl00*abs(v(4)-v(3)); 

v(3)=v2(3); 

v(4)=v2(4); 

end 
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if abs(true_sti£fness(csel(mid__index),cset(mid_index)) - v(3))< .25*abs(v(4)-v(3)) 
v2(3)=v(3)-.5*abs(v(4)-v(3)) ; 
v2(4)=v(4); 

elseifabs(true_stiffness(csel(inid_mdex),cset(mid_iiidcx)) - v(4)) < .25*abs(v(4)-v(3)) 
v2(4)=v(4)+.5*abs(v(4>v(3)) ; 
v2(3)=v(3); 
else 

v2(3)=v(3); 

v2(4>=v(4); 

end 

else 

if abs(loglO((true_stiffness(cset(mid_index),cset(inid_index)) == 0)+... 
tnie_stiffiiess(cset(mid_iiidex),cset(niid_iiidex))) - v(3))< .25*abs(v(4)-v(3)) 
v2(3)=v(3>.5*abs(v(4)-v(3)) ; 
v2(4)=v(4); 

elseif abs(IoglO((tnie_stifIiiess(cset(inid_index),cset(mid_index)) == 0 ) 4 -... 
true_stifiiiess(cset(niid_iudex),cset(mid_index))) - v(4)) < .25*abs(v(4>v(3)) 
v2(4)=v(4)+.5*abs(v(4)-v(3)) ; 
v2(3)=v(3); 
else 

v2(3)=v(3); 

v2(4)=v(4); 

end 

end 

axis(v2) 
grid on 

Iiii=slegend(’COMPUEED STIFFNESS ERROR',TRUE STIFFNESS ERROR'); 
axes(lili) 
if pswitch=='y' 
prtfig(figniim) 
delete(h) 
end 

clear DK 
load DM 
li=tigure(h-H) ; 
fignum=fignuni-Hl ; 
if complete 

plot(\v( 1 :nuiupoiiits-2)/(2*pi),DM(ndx3d([ 1 cset_size*(cset_size+ 1 )/2 nmnpomts-2], 1 ,skycset(niid_index,inid_index),':')),'r— 
»••• 

\v(l:numpoints-2)/(2*pi),tnie_mass(cset(imd_index),cset(niidJndex)).*ones(l,numpoints-2),'b-.') 
ylabeI([T)M(’,int2str(cset(mid_index)),',’,int2str(cset(mid_index)),') lbf-sec''2/in']) 
xlabel(’Omega (Hz)’) 
else 

plot( w( 1 :nuinpoints-2)/(2*pi),Iog 1 0(abs(DM(ndx3d(( 1 csel_size*(cset_size4- 1 )/2 numpoints- 
2], 1 ,skycset(niid_index,iuid_index),':’)))),’r— 

\v( 1 :nuiupoints-2 )/(2 *pi ),log 1 0( . . . 

(... 

true_niass(cset(niid_index),cset(iuid_index)) == 0 ... 

>4-lrue_mass(cset(niid_index),cset(inid_index)). . . 

).*ones( 1 ,numpoints-2),’b-.’) 

yIabel([DM(’,int2str(cset(mid_index)),’,',ml2str(cset(mid_index)),') lbf-sec^2/in (LoglO oO']) 
xlabel(’Omega (Hz)') 
end 
v=a\is; 
if titles 

title('Computed Mass Error vs True Mass Error') 
end 

hold on 

v2(l)=v(l); 

v2(2)=v(2); 
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if complete 
ifabs(v(4)-v(3))< 1 
v2(3)=v(3 > 1 0000 ♦abs( v(4>v(3)); 
v2(4)=v(4)+10000*ab^v(4)-v(3)); 
v(3)=v2(3); 
v(4)=v2(4); 
end 

if abs(true_mass(cset(mid_index),cset(inid_mdex)) - v(3)) < .25*abs(v(4)-v(3)) 
v2(3 )=v(3)-.5 ♦abs(v(4>v(3)); 
v2(4)=v(4); 

elseif abs(true_mass(cset(mid_iiidex),cset(inid_mdex)) - v(4)) < .25*abs(v(4)-v(3)) 
v2(4)=v(4>f.5*abs(v(4)-v(3)); 
v2(3)=v(3); 
else 

v2(3)=v(3); 

v2(4)=v(4); 

end 

else 

if abs(loglO((true_mass(cset(niid_index),cset(mid_index)) == 0)+... 
lrue_niass(csel(mid_index),cset(mid_index))) - v(3)) < .25*abs(v(4)-v(3)) 
v2(3)=v(3>.5*abs(v(4)-v(3)) ; 
v2(4)=v(4); 

elseif abs(logl 0((true__mass(cset(niid_mdex),cset(inid_mdex)) = 0)+. .. 
lnie_mass(cset(inid_mdex),cset(mid_index))) - v(4)) < .25*abs(v(4)-v(3)) 
v2(4)=v(4)^-.5*abs(v(4)-v(3)) ; 
v2(3)=v(3); 
else 

v2(4)=v(4); 

v2(3)=v(3); 

end 

end 

axis(v2) 
grid on 

lili=slegendCCOMPUTED MASS ERROR’,TRUE MASS ERROR’); 
axes(lili) ; 

ifpswitcli=y 
prtfig(fignum) 
delete(h) 
end 

clear DM 

load DC 
li=Qgure(h+l); 
llgnum=fignum+ 1 ; 
if complete 

plot(w( 1 :numpomts-2)/(2*pi),imag(DC(ndx3d([ 1 cset_size*(cset_size+ 1 )/2 numpoints- 
2], 1 ,skycset(inid_index 4 iiid_index),’:’))),'r— 

w( 1 :nuinpoints-2)/(2*pi),iinag(true_damping(cset(mid_index),cset(inid_index)).*ones( 1 4 iumpoiiits- 2 )),'b-.') 
ylabel(['DC(',int2str(cset(mid_index)),V,int2str(cset(niid_mdex)),') Ibf-sec/in']) 
xlabel(’Omega (Hz)’) 
else 

plot(w( 1 mumpoints-2)/(2*pi),logl0(abs(DC(ndx3d([ 1 cset_size*(cset_size+l )/2 nunipoiats- 
2], 1 ,skycset(mid_index,mid_index),’:')))),’r-',.. . 
w( 1 :numpoints-2)/(2*pi),logl 0(abs(. . . 

(... 

true_damping(cset(niid_index),cset(mid_index)) — 0 ... 

)+tnie_dainping(cset(mid_index),cset(mid_index)). . . 

). *ones( 1 4 iumpoints- 2 )),’b-.’) 

ylabel(['DC(',int2str(cset(mid_index)),’,’,int2str(cset(inid_index)),’) Ibf-sec/in (LoglO of)']) 
xlabel(’Oniega (Hz)’) 
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end 

v=axis; 
if titles 

title(’Computed Damping Error vs True Damping Error’) 
end 

v2(l)=v(l); 
v2(2)=v(2); 
if complete 
if abs(v(4)-v(3)) < 1 
v2(3)=v(3)-100*abs(v(4)-v(3)); 
v2(4)=v(4>^-l 00*abs(v(4)-v(3)); 
v(3)=v2(3); 
v(4)=v2(4); 
end 

if abs(true_daniping(cset(mid_mdex),cset(mid_index)) - v(3)) < .25*abs(v(4)-v(3)) 
v2(3)=v(3>.25*abs(v(4)-v(3)) ; 
v2(4)=v(4); 

elseif abs(true_damping(cset(mid_index),cset(mid_index)) - v(4)) < .25*abs(v(4)-v(3)) 
v2(4)=v(4>^-.25♦abs(v(4)-v(3)) ; 
v2(3>=v(3); 
else 

v2(3)=v(3); 

v2(4)=v(4); 

end 

else 

if abs(logl 0((true_damping(cset(mid_mdex),cset(iuid_index)) == 0)^-... 
true__dainping(cset(mid_index),cset(niid_index))) - v(3)) < .l*abs(v(4)-v(3)) 
v2(3)=(v(3)-.25*abs(v(4>v(3))) ; 
v2(4)=v(4); 

elseif abs(loglO((true__damping(cset(niid_index),cset(mid_index)) == Oyi-... 
lrue_dainping(cset(mid_index),cset(mid_index))) - v(4)) < .l*abs(v(4)-v(3)) 
v2(4)=(v(4>^-.25♦abs(v(4)-v(3))) ; 
v2(3)=v(3); 
else 

v2(3)=v(3); 

v2(4)=v(4); 

end 

end 

axis(v2); 
grid on 

Wi=slegend('COMPUTED DAMPING ERROR’,TRUE DAMPING ERROR’); 
axes(lili); 
if ps\vitch==’y' 
prtfig(fignum) 
delete(h) 
end 

clear DC true_damping true_mass true_stifiiiess 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%% 

load H_EXP 
load CORRHA 

%for mid Jndex= 1 :min(cset_size,multplot) 

h=figure(h+l); 

fignimi=fignum+ 1 ; 

plot(\v/(2*pi),log 1 0(abs(H_EXP(ndx3d([ 1 aset_size*(aset_size+ 1 )/2 
numpoints], 1 ,skyred(cset_rel(inid_index),cset_rel(mid_mdex)),':’)))),’r~’,. . . 

\v/(2*pi),log 1 0(abs(CORRHA(ndx3d([aset_size aset_size numpoints], cset_rel(niid_index),cset_rel(mid_index),’:’)))),'b-.') 
ylabel([’H(',int2slr(cset(mid_index)),’,’,int2str(cset(mid_index)),') in/lbf (LoglO of)']) 
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xlabel(’Oniega (Hz)') 
if titles 

titleCExpemnental FRF vs DZ Corrected FRF) 
end 

grid on 
v=axis; 

Wi=slegendCEXPERBvIENTAL FRF/CORRECTED FRF); 
axes(lili) ; 
if psvvitch=='y' 
prtfig(figiumi) 
delete(h) 
end 

clear CORRHA 
load CORRHAD 
h=ligure(h+l); 
figmini=fignuin+l ; 

plot(\v/(2'^pi),log 1 0(abs(H_EXP(ndx3d([ 1 aset_size*(aset_size+ 1 )/2 
luinipoints], 1 ,skyred(cset_rel(mid_index),cset_rel(mid_index)),':')))),'r— . . . 

vv(l :nuiiipoints-2)/(2*pi),loglO(abs(CORRHAD(ndx3d([aset_size aset_size (nuiupoiiits- 
2)],cset__rel(mid_mdex),cset__rel(mid_index),':')))),'b-.') 

ylabel(fH(’,int2str(cset(midJndex)),'/,int2str(cset(inid_index)),') in/lbf (LoglO oO’]) 
xlabel('Oinega (Hz)') 
if titles 

title('Experiniental FRF vs DK/DM/DC Corrected FRF') 
end 
grid on 
v=axis; 

Wi=slegendCEXPERIMElSrTAL FRF','CORRECTED ANALYTIC FRF'); 
axes(lili); 
if pswitch='y' 
prtfig(fignuin) 
delete(h) 
end 

clear CORRHAD 

save plotnum figiumi 
clear 
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%%%%SSTCONF FILE FOR A SPATIALLY COMPLETE BEAM 

beamnidl ;%beamiiidl is a cantilevered 20 element beam model 

%sinmidl ;%simmdl is a masses and springs model 

pswitch=TV ;%%do we print? 

titles=l ;%%display titles 

meters=l ;%use progress meters 

whitebg(lDlack’) ;%switch to black figure background 

close 

lowmode=l; 

liigliniode=10; 

static= 1 ;%%%%%%%%%%%%%reduction metliod 0=Guyan, 1 =IRS,2=Extraction 
%%define A set & O set%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
aset=l :numdof; 



1 5 %oset=2:2:numdof; %%A SPATIALLY INCOMPLETE BEAM%%%%%%%%%% 



%osct=sort(oset); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

oset=[]; %%%%A SPATIALLY COMPLETE BEAM%%%%%%%%%%%%%%% 

20 save sstconf 
sst 

plotsst 
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%%%%SSTCONF FILE FOR A SPATIALLY INCOMPLETE BEAM 

beamindl ;%beamindl is a cantilevered 20 element beam model 

%siinmdl ;%simmdl is a masses and springs model 

pswitch=’n’ ;%%do we print? 

titles=l ;%%display titles 

meters=l ;%use progress meters 

wliitebg(*black') ;%switch to black figure background 

close 

lowmode=l; 

liiglimode=10; 



slatic= 1 ;%%%%%%%%%%%%%reduction metliod 0=Guyan, 1 =IRS,2=Extraction 
%%define A set & O set%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
aset=l rnmndof; 

oset=2:2:numdof; %%A SPATIALLY INCOMPLETE BEAM%%%%%%%%%% 
oset=sort(oset); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



save sstconf 
sst 

plotsst 

beammdl ;%beanmidl is a cantilevered 20 element beam model 

%simmdl ;%simmdl is a masses and springs model 

pswitch=’n' ;%%do we print? 

titles=l ;%%display titles 

meters=l ;%use progress meters 

vvliitebgCblack') ;%switch to black figure background 

close 

lowmode=l; 

liiglunode=10; 

%%define A set & O set%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
aset=l:numdof; 

oset=2:2:nmndof; %%A SPATIALLY INCOMPLETE BEAM%%%%%%%%%% 
oset=sort(oset); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



static=2; 
save sstconf 
dofig4_3 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%% 

%%% COMPARE SINGLE MODE MATRIX AND INTEGRAL %%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%% SOLUTIONS UNDER VARIOUS CONDITIONS %%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
mo%%%%%%%%%% 

clear 

closeall 

loadINT 

whitebg(\vhite’) 

clear L Z_ANAL c_anal c_exp comi h_exp h_anal_red k anal k_exp 
clear kexta kexto ni_anal in_exp mexta mextO lemp_l_diags 
clear lrue_damping true_mass true_stiflhess z_anal z_anal__rcd z_exp 
clear L Z_ANAL c_anal c_exp comi h_exp li_aiial_red k_anal 
clear kexta kexto ni_aiial mexta mextO tenip_l_diags 

loadH.EXP 
load exp 

cset=[5 7 9 11 13 15] 

cset_rel=[3 4 5 6 7 8] 

cset_size=lengtli(cset) 

mid_index=round(lengtli(cset)/2 ); 

flneness=200 

ODZl=[]; 

templ=[]; 

odfreq=omegax( 1 ); 

center_freq=odfreq/(2 * pi) 

wl=odIreq; 

odlengtli=2*pi; 

odivisioiis=2; 

lowerfreq=odfreq-l 5*odleiigtli 
upperfreq=odfreq+25*odlengtli 

owl =odfjreq-,5 *odlength:odleiigtli/odivisioiis:odlieq+. 5 *odlenglli; 
partition=o w 1 /(2 * pi ) 
df=(odlengtli/odi visioiis)/(2 * pi ) 
for count= 1 :lengtli(owl ) 
z_anaCred=kstat-ow 1 (count)^2 * mstat; 
h_aiial_red=iiiv(z_anal_red); 

z__exp=k_exp+j *owl (count) *c_exfK)wl (couiit)^2 *ni_exp; 

h_exp=inv(z_exp); 

h_exp=h_exp(aset,aset); 

hacc=h_anal_red(cset_rel,cset_rel); 

lixcc=h_exp(cset_reLcset_rel); 

dz=inv( inv( inv( hacc)*(hacc-lixcc)* inv(hacc) )-hacc); 

ODZl=[ODZl;dz]; 

temp 1= [tempi; eye(cset_size) -owl(count)^2*eye(cset_size) +j*owl(couiit)*eye(cset_size)]; 
end 

%tlielefldz=ODZl 
%tliebigone=temp 1 
DKDMDCl =templ\ODZl ; 

ODK 1 =DKDMDC 1 ( 1 :cset_size,: ) 

ODMl=DKDMDCl(cset_sizef 1 :2*cset_size,:) 

ODC 1 =DKDMDC l(2*cset_size+l :3 *cset_size,:) 
saveODMl ODMl DKDMDCl ODZl owl wl 
saveODKl ODKl 
save ODCl ODCl 

%i i I I I I l - f H - } M n i M M M M l- f- H -i - f M H t M I I I I t t i M I i | i - t 4 .j.+.t - f i I H M I M M M ) I 

load stat; 
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odcorrha^[]; 

w=lowerfreq:(upperfreq-lowerfreq)/fuieness:upperfreq; 

nuinpoints=lengtli(w); 

kcoirected=kstat; 

nicorrected=nistat; 

ccorrected=cstat; 

kcorrected(cset_rel,cset_rel)=kcorrected(cset_rel,cset__rel)+ODK 1 ; 
mcorrected(cset_rel,cset_rel)=mcorrected(cset_rel,cset_reI)+ODM 1 ; 
ccorrected(cset_rel,cset_rel)=ccorrected(cset_rel,cset_rel)+ODC 1 ; 
for i=l:nuinpoiiils 

teiiipza=kcoiTected+j * w(i )* ccorrected-w(i )^2 *nicorrected; 
tempha=mv(tenipza); 
odcorrha=[odcorrha tempha(:)]; 
end 

imcorrha=[]; 
for i=l:nunipoints 

tempza=kstat+j ^ w(i)*cstat“W(i)^2 *instat; 
teinpha=inv(tempza); 
uncorrha=[unconlia tempha(:)]; 
end 

if meters 

waitbarJiandle=waitbar(0,’Coinputmg Experimental FRF’); 
else 

disp(’Getting experimental FRF) 
end 

H_EXP=[]; 
for count=l:length(w) 
if meters 

waitbai(coimt/length(w)); 

end 

z_exp=k_ex{H-j^w(count)*c_exp-w(count)^2^m_exp; 
h_exp=iiiv(2_exp); 
h_exp=h_exp(aset,aset); 
skyiiidex=l; 
for iiidexl=l:aset_size 
for iiidex2=indexl:aset_size 
red_holder(skyindex)=h_exp(index 1 ,index2 ); 
skyindex=skyiiidex+ 1 ; 
end 
end 

H_EXP=[H^EXP red^holder]; 
end 

if meters 

close(waitbar_handle) 

end 

clear tempha tempza kcoirected mcorrected ccorrected 
figure(l); 

plot(w/(2*pi),loglO(abs(H_EXP(ndx3d([l aset_size*(aset_size+l )/2 nuinpoints],... 
l,skyred(cset_rel(mid_index),cset_rel(inid__index)),’;’)))),’r-,... 

w/ (2 *pi ),log 1 0(abs(odcoiTha(ndx3d( [aset_size aset_size numpoints] ,cset_rel(mid_index),cset_rel(mid_index),':’)))),’g— . . . 
w/(2*pi),logl0(abs(uncorrlia(ndx3d([aset_sizeaset_size numpoints],cset_rel(mid_mdex),cset_rel(mid_index),’:’)))),'b:’) 
ylabel([’H(’,iiit2str(cset(mid_index)),',’,int2str(cset(niid_index)),') iii/lbf (log 10 of)’]) 
xlabel(’Omega (Hz)’) 
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%title('Smgle Mode Coirected FRF (Mode 1 ) ') 
grid 

lili=slegend(’Experimental FRF, Mode 1 Corrected FRF', 'Uncorrectcd Analytic FRF'), 
axes(lili) 

print -dinfile fig5_l 



into dm, dk, & dc %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%% USING MATRDC TECHNIQUE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%Vo%%%%% 
Vo%%%%%%%%Vo%%%%%%%%%%%Vo%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Vo%%%Vo%%% 

tenipl=[]; 

odlengtli=2*pi; 

owl=omegax(l )-5*odlengtli:(omegax(2)-omegax( 1 )+10*odlengtli)/14:omegax(2)+5*odlengtli; 

for count= 1 :lengtli(owl ) 
z_anal_red=kstat-owl(count)^2*mslat; 
h_anal_red= in v(z__anal_red); 

z_exp=k_exp+j *owl (coimt)*c_exp-owl (count)^2 *m_exp; 

h_exp=inv(z_exp); 

h__exp=h__exp(aset,aset); 

liacc=h_anal_red(cset__rel,cset__rel); 

lixcc=h__exp(cset_rel,cset_rel); 

dz=inv(inv(inv(hacc)*(liacc-hxcc)*iiiv(hacc))-hacc); 

ODZl=[ODZl; dz]; 

temp 1= [tempi; eye(cset_size) -owl(count)'^2*eye(cset_size) +j*o\vl (count)* eye(cset_size)]; 
end 

DKDMDCl=tenipl\ODZl ; 

ODKl =DKDMDC 1 ( 1 xset^size,:) ; 

ODMl=DKDMDCl(cset_size^-l:2*cset_size,:) ; 

ODCl =DKDMDC 1 (2*cset_size^-l :3*cset_size,:); 
save ODMl ODMl DKDMDCl ODZl owl wl 
save ODKl ODKl 
save ODCl ODCl 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

imm^/o%%%%Vo%% 

odcoirha=[]; 

w=omegax(l)-l 5*pi:(omegax(2)-omegax( 1 )+35 *pi)/fineness:oinegax(2)+20*pi; 

numpoints=lengtli(w); 

kcorrected=kstat; 

mcorrected=mstat; 

cconrected=cstat; 

kcorrected(cset__rel,cset_rel)=kcoiTected(cset__rel,cset_rel)+ODK 1 ; 
mcoirected(cset__reL,cset_rel )=mcoirected(cset_rel,cset_rel)+ODM 1 ; 
ccorrected(cset_rel,cset_rel)=ccoiTected(cset_reL,cset_rel)+ODC 1 ; 
for i=l:numpoints 

tempza=kcorrected+j * w(i)*ccorrected-w(i)'^2 *nicorrected; 
lempha=inv(tempza); 
odcorrlia=[odcorrha tempha(:)]; 
end 

micorrha=[]; 
for i=l:numpoints 

teinpza=kstat+j * w(i)*cstat-w(i)^2 *mstat; 
tempha=inv(terapza): 
uncorrha=[uncorrha tempha( : )] , 
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end 

if meters 

\vaitbarJiaiidle=waitbar(0,'Computiiig Experimental FRP')i 
else 

disp('Getting experimental FRF) 
end 

H^EXP=[]; 
for count= 1 ;lengtli(w) 
if meters 

waitbar(count/lengtli(w)); 

end 

z_exp=k_exp+j * \v(coimt)* c_exp-\v( coimt)^2 * m_exp; 
h_exp=inv(z_exp); 
h_exp=h_exp(aset,aset); 
skyindex=l; 
for indexl=l:aset_size 
for mdex2=indexl:aset_size 
red Jiolder(slcyindex)=h_exp( index 1 ,index2); 
skyindex=skyindex+ 1 ; 
end 
end 

H_EXP=[H_EXP red^^holder^]; 
end 

if meters 

close( waitbar_handle) 
end 

clear tempha tempza kcorrected mcorrected ccoixected 
figure(2); 

plot(\v/(2*pi),loglO(abs(H_EXP(ndx3d([l aset_size*(aset_size+l )/2 numpoints],... 

1 ,skyred(cset_rel(mid_index),cset_rel(mid_mdex))/:’)))),’r-’,. . . 

\v/(2*pi),logl0(abs(odcoirha(ndx3d([aset_size aset_size numpoints],cset_rel(mid_index),cset_rel(mid_index),V)))),*g-',... 
\v/(2*pi),loglO(abs(uncorTha(ndx3d([aset__sLze aset_size numpoints], cset_rel(mid_index),cset_rel(mid_index),*:’)))),l):’) 
ylabel(['H(*,int2str(cset(inid_index)),’,’,int2str(cset(mid_index)),') in/lbf (loglO oQ']) 
xlabel(’Omega (Hz)') 

%title('Single Mode Corrected FRF (Mode 1 ) *) 
grid 

lili=slegend(Experimental FRF’, 'Single Mode Corrected FRF, 'Uncorrected Analytic FRF); 
axes(lili) 

print -dm file fig5_2 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
into ODM, ODK, & ODC %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
^t^/li5^%Jsft64^4i^^C&eq-lo\verfreq)/fineness:up^^ 
lenctrl=[25 10 1]; 
odivisions=3; 
for ind= 1 :lengtli(lenctrl) 
odlengtli=lenctrl(ind)*2*pi; 

ODZ=[j; 

temp=[]; 

ow=odfreq-.5 *odleiigtli:odlengtli/odi visions:odIreq+.5 *odlength; 
for comit=l :lengUi(ow) 
z_anal_red=kstat*ow( count )^^2 * mstat; 
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h_aiial_red=iii v(z_anal_red ); 

z__exp=k_exp+j *ow(count)*c_exp-ovv(coimt)^2 *m_exp; 

h_exp=inv(z_exp); 

h_exp=h_exp(aset,aset); 

hacc=h_anal_red(cset_rel,cset_rel); 

lixcc=h_exp(cset_rel,cset_rel); 

dz=mv(inv(inv(hacc)^(hacc-lLxcc)*inv(hacc))-hacc); 

ODZ=[ODZ; dz]; 

teinp=[temp; eye(cset_size) -ow(count)^2*eye(cset_size) +j*ow(count)*eye(cset_size)]; 
end 

DKDMDC=temp\ODZ; 

ODK=DKDMDC(l:cset_size,:) ; 

ODM=DKDMDC(cset_size+'l:2*cset_size,:) ; 

ODC=DKDMDC(2*cset_size+l :3*cset_size,:); 

clear hacc lixcc h_exp z_exp z_anal_red h_aiial_red dz 

save ODM ODM DI^MDC ODZ ow 

save ODK ODK 

save ODC ODC 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

waitbar_handle=waitbar(0,’Computing Experimental FRF); 
else 

disp(’Getting experimental FRF) 
end 

H^EXP=G; 

for count=l :lengtli(plotvv) 
if meters 

waitbar(count/lengtli(plotw)); 

end 

z_exp=k_exp+-j*plotw(count)*c_exp-plotw(couiit)^2*m_exp; 
h_exp=inv(z_exp); 
h_exp=h_exp(aset,aset); 
skyindex=l; 
for index 1 = 1 :aset_size 
for index2=indexl :aset_size 
red_holder(skyindex)=h_exp( index 1 ,index2 ); 
skyindex=skyindex+ 1 ; 
end 
end 

H_EXP=[H_EXP redjiolder’]; 
end 

if meters 

close(waitbar_handle) 

end 

odcorrha=[]; 

kcorrected=kstat; 

mcorrected=mstat; 

ccorrected=cstat; 

kcorrected(cset_rel ,cset_rel )=kcorrected(cset_rel ,cset_rel)+ODK; 
mcorrected(cset_rel,cset_rel)=mcorrected(csel_rel,cset_rel)+ODM; 
ccorrected(cset_rel,cset_rel)=ccorrected(cset_rel,cset_rel)+ODC; 
for i=l:lengtli(plot\v) 

tempza=kcorrected+j*plotvv(i)*ccorrected-plotw(i)^2*mcorrected; 
tempha=inv(tempza); 
odcorrlia=[odcorrha tempha(:)]; 
end 

if ind=l 
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plotl=odcorrha; 
elseif ind==2 
plot2=odcorrha; 
elseif iiid—3 
plot3=odcorrha; 
else 

plot4=odcorrha; 

end 

end 

figure(3); 

plot(plotw/(2*pi),logl0(abs(H_EXP(iidx3d{[l aset_si2e*(aset_size+l)/2 lengtli(plotw)],... 

1 ,skyred(cset_rel{mid_index),cset_rel(mid_index)),':')))),'r- . . 

plotw/(2*pi), Iogl0(abs(plotl(ndx3d([aset_si2e asel_size leiigtli(plotw)],cset_rel(inid_index),cset_rel(mid_index),’:'))))/g— 

plotw/(2*pi), Iogl0(abs(plot2(ndx3d([aset_si2e aset_size 
lengtli{plotw)],cset_rel(inid_mdex),cset_rel(niid_index),V)))),T?:',... 

plotw/(2*pi), log 1 0(abs(plot3(ndx3d([aset_size aset_size lengtli(plotw)],cset_rel(inid_index),cset_rel(mid_iiidex),':')))),’k-.') 

ylabel(fH(’,int2str(cset(imd_index)),',',iiit2str(cset(mid_iiidex)),’) in/lbf (loglO of)']) 

xlabel('Omega (Hz)') 
grid on 

%titleCEXPERMENTAL FRF AND CORRECTED FRFS USING MODE 5 SOLUTIONS’) 
lili=legend(Exp FRF', ... 

[nmu2str(lenclrl(l)),’ Hz Bandwidtli Corrected FRF'], ... 

[nuni2str(lenctrl(2)),' Hz Bandwidtli Corrected FRF'], ... 

[num2str(lenclrl(3)),' Hz Bandwidtli Corrected FRF]) ; 
axes(hli) 

print -diiifile fig5_3 



odlengtli=l*2*pi; 
countctrl=[3 10 50 200] ; 
for ind=l:lengtli(countctrl) 
odi visions=countctrl(iiid> 1; 

ODZ=[]; 

tenip=[]; 

ow=odfreq-. 5 ♦odlengtli:odlengtli/odi visions:odfreq+. 5 ^odlengdi; 
for count=T.lengtli(ow) 
z_anal_red=kstat-ow(count)'^2 *iiistat; 
li_aiial_red=iiiv(z_anal_red); 

z_exp=k_exp+j *o w(count)*c_exp-o w(count)'^2 *m_exp; 

li_exp=inv(z_exp); 

li_exp=li_exp(aset,aset); 

hacc=li_anal_red(cset_rel,cset_rel); 

lixcc=li_exp(cset_rel,cset_rel); 

dz=inv(inv(inv(hacc)*(liacc-lL\cc)*inv(liacc))-hacc); 

ODZ=[ODZ; dz]; 

temp=[temp; eye(cset_size) -ow(count)^2*eye(cset_size) +j*ow(comit)*eye(cset_size)]; 
end 

DKDMDC=tenip\ODZ; 

ODK=DKDMDC(l:cset_size,:) ; 

ODM=DKDME>C(cset__size+ 1 :2 *cset_size,: ) ; 

ODC=DKDMDC(2*cset_size+l :3 *cset_size,:); 
clear hacc lixcc li_exp z_exp z_anal_red li_anal_red dz 
save ODM ODM DKDMDC ODZ ow 
saveODKODK 
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save ODC ODC 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%% 

odcorrha=[]; 

kcorrected=kstat; 

mcx)iTected=mstat; 

ccorrected=cstat; 

kcorrected(cset_rel,cset_rel)=kcorrected(cset_rel,cset_rel)+ODK; 
mcorrected(cset_rel,cset_rel)=nicx)iTected(cset_reKcset_rcl>i-ODM; 
ccorrected(cset_rel,cset_rel)=ccoiTected(cset_rel,cset_rel)+ODC; 
for i=l :lengtli(plotw) 

tenipza=kcorrected+j *plotw( i )*ccorrected“plotvv(i *incorrected; 

tenipha=iii v( teinpza); 
odcorrha=[odcorrha teinpha(:)]; 
end 

if ind==l 
plotl=odcorrha; 
elseif ind=2 
plot2=odcorrha; 
elseif ind=3 
plot3=odcoirha; 
else 

plot4=odcorrha; 

end 

end 

figure{4); 

plot(plotw/(2*pi),logl0(abs(H_EXP(ndx3d([l aset_size*(aset_size+l )/2 lengtli(plotw)],... 
l,skyred(cset_rel(midjndex),cset_rel(midjndex)),':')))),'r-,... 

plot\v/(2*pi), Iogl0(abs(plotl(ndx3d([aset_si2e aset_size lengtli(plotw)],cset_iel(niid_index),cset_rel(mid_index),’:')))),'g- 

plotw/(2*pi), Iogl0(abs(plot2(ndx3d([aset_size aset_size 
lengtli(plotw)],cset_rel(inid_index),cset_rel(mid_iiidex),’:’)))),'b:',... 

ploUv/(2*pi), Iogl0(abs(plot3(ndx3d([aset_si2e aset_size length(plotw)],cset_rel(niid_index),cset_rel(mid_index),’:’)))),’k- 

plot\v/(2*pi), Iogl0(abs(plot4(ndx3d([aset_size aset_size length(plotw)],cset_rel(niid_index),cset__rel(inidJndex),':’)))),'iuO 

ylabel(('H(’,int2str(cset(inid_index)),V,int2str(cset(inid_index)),’) in/lbf (loglO of)']) 

xlabel(’Omega (Hz)') 
grid on 

%titleCEXPERMENTAL FRF AND CORRECTED FRFS USING MODE 1 SOLUTIONS') 
lili=legend(’Exp FRF’, ... 

(nuin2str(countctrl(l)),' Pts Corrected FRF'], ... 

[nuni2str(countctrl(2)),' Pts Corrected FRF], ... 

(nmn2str(countctrl(3)),' Pts Corrected FRF'], ... 

(nuni2str(countctrl(4)),’ Pts Corrected FRF]); 
axes(lili) 

print -dinfile fig5__4 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
into dm, dk, & dc %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%% USING INTEGRAL TECHNIQUE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%plot\v=lowerfreq:d\v:upperfreq; 

lcnctrl=[25 10 1]; 

odivisions=3; 

for ind= 1 :lengtli(lenctrl) 
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odlengtli=lenctrl(iiid)*2^pi; 

owl =odfreq-. 5 *odlengtli:odleiigtli/(odivisions);odfreq+. 5 ^odlengtli; 

Wl=owl; 

%%coinpute integrals %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%weiglit=ones(size(owl ))*(omegax( 1 )); 

weiglit=owl; 

DZ^R=[]; 

DZ_I=[]; 

for coimt=l :lengtli(owl ) 

z_anal_red=kstat/(j *owl (coimt))+j *ow 1 (count)*mstat; 
h_aiial_red=mv(z_aiial_red); 

z_exp=k_exp/0‘ *owl (count))+c_exp/(j *owl (count ))+j *owl (count) *m_exp; 

h_exp=inv(z_exp); 

h_exp=h_exp(aset,aset); 

hacc=h_anal_red(cset_rel,cset_rel); 

lixcc=h_exp(cset_rel ,cset_rel); 

dz=inv(inv(inv(liacc)*(hacc-hxcc)*inv(hacc))-hacc); 

dz_r=real(dz); 

dz__i=imag(dz); 

DZ_R=pZ_Rdz_r(:)]; 

DZ_I=PZJ dzj(:)]; 
end 

[INTK INTM INTC]=intsub(DZ_I,DZ_R,owl ,W 1 , 1 ./(j *owl ),cset_size); 

odcorrha=[]; 

load stat 

kcorrected=kstat; 

mcorrected=mstat; 

ccorrected=cstat; 

kcorrected(cset_rel,csel_rel)=kcorrected(cset_rel,cset_rel)+INTK; 
incorrected(cset_rel,cset_rel)=incoiTected(cset_rel,csel_rel)+INTM; 
ccorrected(cset_rel,cset_rel)=ccorrected(cset_rel,cset_rel)+INTC; 
for i=l :lengtli(plotw) 

tempza=kcorrected+j*plotw(i)*ccoiTected-plotw(i)^2*mcoiTected; 
tempha=inv(tempza); 
odcorrha=[odcoiTha tempha(:)]; 
end 

if ind== 1 

intplotl=odcorrha; 
save INTMl INTM 
save INTKl INTK 
save DSTTCl ENTC 
elseif ind==2 
intplot2=odcorrha; 
save INTM2 INTM 
save INTK2 INTK 
save INTC2 INTC 
elseif ind=3 
intploG=odcorrha; 
save INTM3 INTM 
save INTK3 INTK 
save INTC3 INTC 
else 

intplot4 =odcorrha; 
save INTM4 INTM 
save INTK4 INTK 
save ENTC4 INTC 
end 
end 
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if meters ' 

waitbar_handle=waitbar(0,'Computiiig Experimental FRF); 
else 

dispCGetting experimental FRF’) 
end 

H_EXP=[]; 

for coimt=l:lengtla(plotw) 
if meters 

waitbar(count/lengtli(plotw)); 

end 

z_exp=k_exp+j*plotw(coimt)*c_exp-plotvv(count)^2'^m_exp; 
h__exp=inv(2_exp); 
h_exp=h_e\p(aset,aset); 
skyindex=l; 
for indexl = 1 :aset_size 
for index2=mdexl:aset_si2e 
red_holder(skyindex)=h_exp(mdex 1 ,index2); 
skymde\’=skyindex+ 1 ; 
end 
end 

H_EXP=[H_EXP redjioldei"]; 
end 

if meters 

close(waitbar_haiidle) 

end 

figure(5); 

plot(plotw/(2*pi),loglO(abs(H_EXP(ndx3d([l aset_size*(aset_size+l )/2 length(plotw)],... 

l,skyred(cset_rel(mid_index),cset_rel(mid_index)),V)))),'r-V- 

plotw/(2'*‘pi), Iogl0(abs(intplotl(ndx3d([aset_size aset_size 
lengtli(plotw)] ,cset_rel(mid_index),cset_rel(mid_index),V)))),’g- . . 

plotw/(2*pi), Iogl0(abs(mtplot2(ndx3d([aset_si2e aset_size 
lengtli(plotw)],cset_rel(mid_index),cset_rel(mid_index),':')))),’b:',... 

plotw/(2*pi), Iogl0(abs(mtplot3(ndx3d([aset_si2e aset_size 
length(plotw)],cset_rel(mid_index),cset_rel(mid_index),V)))),'k-.') 
ylabel(rH(',int2str(cset(mid_mdex)),V4i^f2str(cset(mid_iiidex)),’) in/lbf (log 10 of)’]) 

xlabel(’Omega (Hz)') 
grid on 

%titleCEXPERMENTAL FRF AND CORRECTED FRFS USING MODE 1 INT SOLUTIONS') 
liIi=legendCExp FRF', ... 

[nnm2str(lenctrl(l)),' Hz Bandwidtli Corrected FRF'],... 

[num2str(Ienctrl(2)),' Hz Bandwidtli Corrected FRF],... 

[num2str(lenctrl(3)),' Hz BandwidUi Corrected FRF']); 
axes(lili) 

print -dmfile fig5_5 

odlengtli=2*pi; 
countctrl=[3 10 50 200]; 
for ind=l:lengtli(coimtctrl) 
odivisions=countctrl(ind); 

owl =odfi*eq-. 5 * odlengtli:odlengtli/(odi visions-1 ):odfreq+. 5 *odlengtli; 

Wl=owl; 

%%compute integrals %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%weiglit=oiies(size(owl ))*(omegax( 1 )); 
weiglit=owl; 
if ind == 1 
owl 

df=(odlengtli/(odi visions- 1 ))/(2 ^ pi) 
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end 

%% compute integrals %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%o/o%%%%o/oO/^,o/^ 
DZ_R=Q; 

DZ_I=Q; 

for count= 1 :lengtli(owl ) 

z_anal_red=kstat/(j *o w 1 (couiit)>4'j *owl (count)*mstat; 
h_anal_red=inv(z_aiial_red); 

z_exp=k_exp/(j *o wl (count))+c_exp/(j *owl (count))+j *owl (couiit)*m_exp; 

h_exp=inv(z_exp); 

h_exp=h_exp(aset,aset); 

hacc=h_aiial_red(cset_rel,cset_rel); 

hxcc=h_exp(cset_rel,cset_rel); 

dz=inv(inv(mv(hacc)*(hacc-lixcc)*inv(hacc))-hacc); 

dz_r=real(dz); 

dzj=iinag(dz); 

DZ_R=PZ_R dz_r(:)]; 

DZJ=PZ_I dz_i(:)]; 
end 

if ind == 1 
DZ_R; 

DZ^t 

for i=l:lengtli(owl) 

dztempl=DZ_R(ndx3d([ 1 cset_size*(cset_size+ 1 )/2 lengtli(owl )], 1 ,1 :cset_size*(cset_size+l )/2,i)); 

dzteinp2=DZ_I(ndx3d([l cset_size*(cset_size+l )/2 lengtli(o\vl)],l,l :cset_size*(cset_size+l )/2,i)); 

skyindex=l; 
for index 1 = 1 xset size 
for index2=index 1 :cset_size 
dz I ( index 1 ,index2)=dztemp 1 (sky mdex); 
dz 1 (index2,index I )=dztempl (skyindex); 
dz2(index 1 ,mdex2)=dztemp2(skyindex); 
dz2(index2,index I )=dzteinp2( skyindex); 
skyindex=skyindex+l ; 
end 
end 
dzl 
dz2 
end 

end 

[USTTK INTM INTC]=mtsub(DZ_ipZ_R,o\vl,Wl,l./(j*o\vl),cset_size); 
ifind== I 

[INTK INTM INTC]=mdsub(DZ_I,DZ_R,owl,Wl,l./(j*owl),cset_size); 

INTM 

INTK 

INTC 

end 

ydmmmff/0%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

mi^stat 

kcorrected=kstat; 

nicorrected=instat; 

ccorrected=cstat; 

kcorrected(cset_rel,cset_rel)=kcorrected(cset_rel,cset_rel)+INTK; 
mcoirected(cset_rel,cset_rel)=incorrected(cset_rel,cset_rel)4-INTM; 
ccorrected( cset_rel ,cset_rel)=ccorrected(cset_rel ,cset_rel )+INT C ; 
for i= 1 :lengtli(plotw) 

tempza=kcorrected+j *plotw(i )*ccorrected-plotw(i)^2 *incorrected; 
tempha=inv(teinpza); 
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odcorrha= [odcorrha tempha( : )] ; 
end 

if ind==l 

intplotl=odcorrha; 
save INTMl INTM 
save INTKl INTK 
saveINTCl INTC 
elseif iiid==2 
intplot2=odcorrha; 
save INTM2 INTM 
saveINTK2 INTK 
save INTC2 INTC 
elseif ind==3 
intploO =odcorrha ; 
save INTM3 INTM 
save INTK3 INTK 
save INTC3 INTC 
else 

intplot4=odcorrlia; 
save INTM4 INTM 
save INTK4 INTK 
save INTC4 INTC 
end 
end - 

if meters 

\vaitbar_handle=waitbar(0,'Conipiitiiig Experimental FRF); 
else 

disp(’Gettmg experimental FRF') 
end 

H_EXP=[]; 

for coimt=l:lengtli(plotw) 
if meters 

waitbar(count/lengtli(plotw)); 

end 

z_exp=k_exp+j *plotw(count)*c_exp-plotw(count)'^2 *ni_exp; 
h__exp=inv( z_exp); 
h_exp=li_exp(aset,aset); 
skyindex=I; 
for index I = 1 :aset_size 
for index2=indexl :aset_size 
red_holder(skyindex)=h_exp(index 1 ,index2); 
skyindex=skyindex+ 1 ; 
end 
end 

H_EXP=[H^EXP red^holder']; 
end 

if meters 

close( wai tbar_hand le) 
end 

figure(6); 

plot(plotw/(2*pi),logl0(abs(H_EXP(ndx3d([l aset_size*(aset_size+l )/2 lengtli(plotw)],... 

1 ,sk>Ted(cset_rel(mid_index),cset_rel(midJndex)),V)))),'r-',. . . 
plotw/(2*pi), Iogl0(abs(intplotl(ndx3d([aset_size aset_size 
lengtli(plotw)],cset_rel(mid_index),cset_rel(mid_index),’:’)))),’g-',... 

plolw/(2*pi), Iogl0(abs(intplot2(ndx3d((aset_size aset_size 
lengtli(plotw)],cset_rel(mid_index),cset_rel(mid_index),’:')))),'b;',... 
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plotvv/(2*pi), Iogl0(abs(intplot3(ndx3d([aset_size aset_size 
lengtli(plotvv)] ,cset_rel(niid_index),cset_rel(mid_index),':’)))),'k-.',. . . 

plotw/(2*pi), Iogl0(ab$(intplot4(ndx3d([aset_size aset_size 
leiigtli(plot\v)],cset__rel(mid_iudex),cset_rel(niid_mdex),’:’)))),’m:') 
ylabel([’H(’,int2str(cset(mid_mdex)),V,int2str(cset(inid_index)),’) iii/lbf (log 10 of)’]) 

xlabel('Oinega (Hz)') 
grid on 

%tiaeCEXPERMENTAL FRF AND CORRECTED FRFS USING MODE 1 INT SOLUTIONS') 
liii=slegendCExp FRF, ... 

[num2str(coimtctrl(l)),' Pts Corrected FRF'], ... 

[num2str(countctrl(2)),' Pts Corrected FRF'], ... 

[num2str(countctrl(3)),' Pts Corrected FRF], ... 

[num2str(countctrl(4)),’ Pts Corrected FRF]); 
axes(Wi); 

print -dnrfile fig5_6 
figure(7); 

plot(plotw/(2^pi),loglO(abs(H_EXP(ndx3d([l aset_size*(aset_size+l )/2 lengtli(plotvv)],... 

1 ,skyred(cset_rel(inid_mdex),cset_rel(nud_mdex)),':’)))),'r“', . . . 

plotw/(2*pi), Iogl0(abs(plotl(ndx3d([aset_sizeaset_size lengtli(plotw)],cset_rel(mid_index),cset_rel(mid_index),':')))),'g- 

f 

plot\v/(2^pi), Iogl0(abs(intplotl(ndx3d([aset_size aset_size 
lengtli(plotw)],cset_rel(mid_index),cset_rel(inid_index),':')))),'b:') 
ylabel(['H(',int2str(cset(mid_index)),V,uit2str(cset(mid_mdex)),') in/lbf (log 10 of)']) 

xlabel('Omega (FIz)') 
if titles 

% tideCEXPERMENTAL FRF AND CORRECTED FRFS USING MODE 1 OD&INTGRL SOLTNS') 
end 
grid on 

lili=legendCExp FRF, ... 

[nuin2str(lenctrl(l)),' Hz Bandwidtli MAT Corrected FRF'], ... 

[num2str(leiictrl(l)),' Hz Bandwidtli INT Corrected FRF]); 
axes(lili) 
if pswitch=='y’ 
print -dcdjcolor 
end 

print -dinfile fig5_7 
figure(8); 

plot(plotw/(2*pi),loglO(abs(H_EXP(ndx3d([l aset_size*(aset_size+l )/2 length(ploriv)],... 
l,skyred(cset_rel(mid_index),cset_rel(niid_index)),':')))),'r-',... 

plotw/(2 *pi ), log 1 0(abs(plot2(ndx3d([aset_size aset_size lengtli(plotw)] ,cset_rei( inid_index),cset_rel( mid_index),':' )))),'g— 

I 

,... 

plotw/(2*pi), Iogl0(abs(intplot2(ndx3d([aset_size aset_size 
lengtli(plotw)],cset_rel(niid_index),cset_rel(mid_mdex),':’)))),'b:’) 
ylabel( rH(',int2str(cset(mid_index)),’,',int2str(cset(niid_index)),') in/lbf (log 1 0 of)'] ) 

xlabel('Omega (Hz)') 
if titles 

%titleCEXPERMENTAL FRF AND CORRECTED FRFS USING MODE 1 OD&INTGRL SOLTNS') 
end 
grid on 

lili=legend(Exp FRF, ... 

[nuni2str(lenctrl(2)),' Hz Bandwidth MAT Corrected FRF'], ... 

(iium2str(lenctrl(2)),' Hz Bandwidtli INT Corrected FRF']); 
axes(hli) 
if ps\\itch=='y' 
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print -dcdj color 
end 

print "dinfile fig5_8 
figure(9); 

plot(plotw/(2*pi),loglO(abs(H_EXP(n<ix3d([l aset_size*(aset_size+l )/2 lengtli(plotvv)],... 

1 ,skyred(cset_rel(mid_index),cset_rel(niid_index)),’:’)))),'r- . . 

plotw/(2*pi), log 1 0(abs(plot3(ndx3d([aset_size aset_size lengtli(plotvv)],cset_rel(mid_index),cset_rel(mid_index),':’)))),'g- 

plot\v/(2*pi), Iogl0(abs(intplot3(ndx3d([aset_size aset_size 
lengtli(plot\v)] ,cset_rel(inid_index ),cset_rel(niid_index),':’)) )),’b:’) 
ylabel([’H(’,int2str(cset(inid_mdex)),V,iiit2str(cset(mid_mdex)),') in/lbf (loglO oQ’]) 

xlabel(’Omega (Hz)') 
if titles 

%titleCEXPERMENTAL FRF AND CORRECTED FRFS USING MODE 1 OD&DSTTGRL SOLTNS’) 
end 
grid on 

lili=legendCExp FRF’, ... 

[nuni2str(lenctrl(3)),’ Hz Bandwidtli MAT Corrected FRF'], ... 

[nuni2str(lenc^l(3)),' Hz Bandwidtli INT Corrected FRF]); 
axes(hli) 
if pswitch='y' 
print -dcdjcolor 
end 

print -dnifile fig5_9 
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clear . 
odivisions=3 
startloop=l 
skiploop=l 
5 eiidloop=4 

cd incomp 1 

save odconf odivisioiis startloop skiploop endloop 
chap6_l 
chap6_2 
10 cd.. 

clear 

odivisions=l 
startloop=3 
skiploop=l 
1 5 endloop=4 

cd compl 

save odconf odivisioiis startloop skiploop endloop 

chap6_l 

chap6_3 
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%%%CHAP6J%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%% 

%%%Deconipose DZ into DM, DK, & DC%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%using matrix formulation%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%% 
clc 
clear 
closeall 
loadINT 
%if pswitch==y 
wliitebg(’wliite') 
close 
%end 
h=0; 

fignum=l; 

use_aiitires=l; 

litles=0; 

imd_index=round(lengtli(cset)/2 ); 

clear L Z_ANAL c_exp conn h_anal_red 
clear kexto m_exp mexto temp__l_diags 

clear true_damping true_mass true_stiffiiess z_anal z_anal_red z__exp 
clear L Z_ANAL c_exp conn h_anal_red 
clear kexta kexto mextO temp_l_diags 

%%%%%%%%%%%%%%%%%%%FORCETHECSET%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%SPAT1ALLY INCOMPLETE%%m%%%%%%m%%%%%ro 
if complete 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%SPATIALLY COMPLETE%%%%%%%%%%%%%%%%% 
cset=[4 5 6 7 8 9 10 1 1]%%%%%%%%%%%%%%%%%%%%% 
csetj:el=csct%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
cset_sizc=kngth(cset)%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
else 

%%%%%%%%%%%%SPATIALLY INCOMPLETE%%%%%%%%%%%%%%%%%% 
cset=[l 3 5 7 9 11 13 15]%%%%%%%%%%%%%%%%%%%%%% 
cset_rel=[l 2 3 4 5 6 7 8]%%%%%%%%%%%%%%%%%%%%% 
cset_size=lengUi(cset)%%%%%%%%%%%%%%%%%%%%%%%%% 
end 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



load exp 
odivisions=3 ; 
load odconf 
nuinniodes=8 ; 
odlengdi=2.00 1 *pi ; 
firsttime=l; 

dw=odlengtli/(odi visions^- 1 ); 
do\v=(freqtop-freqbottoiu)/( fineness- 1 ); 
owref=freqbottom :do w:freqtop; 

if firsttime 
if meters 

vvaitbar_handle=waitbar(0,’Computing Experimental FRF'); 
else 

dispCGettiiig experimental FRF') 
end 
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H_.T^EXP=U; 
for count=l:lengtli(owref) 
if meters 

waitbar(count/lengtli(owref)); 

end 

z_exp=k_exp+j *owref{count)^c_exfK) wref(couiit)^2 *m_exp; 

h_exp=inv(z_exp); 

h_exp=h_exp(aset,aset); 

% skyiiidex=l; 

% for index 1 = 1 :aset_size 
% for index2=mdex 1 :aset_size 
% red_holder(skyindex)=h_exp(index 1 ,index2 ); 

% skyindex=skyiiidex+ 1 ; 

% end 
% end 

% H^T_EXP=[H_T_EXP redjiolder’]; 

H_T_EXP=[H_T_EXP h_exp(:)]; 
end 

if meters 

close(waitbar_handle) 

end 

saveJH_T^EXP H_T_EXP 
clear H^T^EXP 
%loadH_ANAL_RED; 

%uncorrha=H_ANAL_RED; 

%clear H^ANAL^RED 
uncorrha=[]; 
if meters 

vvaitbar_handle=waitbar(0,'Computing Uncorrected FRF'); 
else 

disp('Getting Uncorrected FRF') 
end 

load stat 

for i=l:lengtli(owref) 
if meters 

waitbar(i/lengtli(o\vref)); 

end 

tempza=kstat+j *owref(i)*cstat-owref(i)^2 *mstat; 
tempha=inv(tempza); 
uncorrha=[uncoiTha tempha(:)]; 
end 

if meters 

close(waitbarJiandle) 

end 

save imcorrha uncoirha 
clear uncorrha 
end 

loadH_T^EXP 

temp=[); 

for index= 1 munimodes 
for coord=l :aset_size 
if index == 1 

antiresfreqs=fmd(owref>0 & owref<oniegax(index)); 
else 

antiresfreqs=fuid(owref>omegax( index- 1 ) & owref<omegax( index)); 
end 

antires=min(log 1 0(abs(H_T_EXP( ndx3d( [aset_size aset_size lengtli(ovvref)] , • • • 
coord,coord,antiresfreqs))))); 

\vlire=rmd(logl 0(abs(H_T_EXP(ndx3d([aset_size aset_size lengtli(o\\TeO],.-- 
coord,coord,antiresfreqs))))==antires); 
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teinp=o\vref(aiitiresfreqs) ; 
antiresfreq(coord,index)=temp(wlire( 1 )); 
end 
end 

temp=[]; 

for index=l:nummodes 

resfreqs=find(owref^omegax(index)-dow & owref<omegax( index )+'dow); 
res=max(log 1 0(abs(H_T__EXP(ndx3d( [aset_size aset_size length(o\vreO] ■ . 
cset_rel( 1 ),cset_rel( 1 ),resfreqs))))); 

\vlire=find(logl 0(abs(H_T_EXP(ndx3d( [aset_size aset_size lenglli(o\vret)],. .. 

cset_rel( 1 ),cset_rel( 1 ),resfreqs))))==res); 
temp=owref(resfreqs) ; 
resfreq(index)=tenip(wlire( 1 )); 
end 

clear H_T__EXP 

%clear omegax 
%omegax=resfreq; 

for loopmdex=startloop:skiploop;endloop 
ODZ=[]; 
temp=[]; 
rdiag=[]; 
ow=[]; 

lo\\'\ver=o\vref( 1 )-.5*odlengtli; 
uppper=owref(l>+'.5*odlengtli ; 

% if odi visions == 1 
% o \v=[ow owref( 1 )] 

% else 

% ow=[ow lowwer+d\v:dw:uppper-dw]; 

% end 

% rdiag=[rdiag ones( 1 ,odivisions*cset_size). *o\vref( 1 )]; 



for index=l:loopindex 
% lowwer=omega\(mdex)-.5*odlengtli ; 

% uppper=omegax(index)+'.5*odlength ; 

lo\v\ver=resfreq(index)-. 5 *odlaigtli ; 

uppper=resireq(index>+'.5*odlengtli ; 

o\v=[ow lowwer+dw:d\v;uppper-dw]; 

rdiag=[rdiag ones(l ,odivisions*cset_size).*resfreq(index)]; 

if use_antires 

for coord=l;aset_size 

lowwer=antiresfreq(coord,index)-. 5 *odlengtli; 
uppper=antiresfreq(coord,index)+-.5*odIengtli; 
size(lowwer+dvv:dw:uppper-dw) 
ow=[ow lowwer+dw:dw:uppper-dw] 
rdiag=[rdiag ones( 1 ,odivisions*cset_size). *resfreq( index)]; 
end 
end 
end 

R=diag(rdiag); 
load stat 
if meters 

\vaitbar_handle=waitbar(0,'Coniputmg DK/DM/DC’); 
else 

dispCGetting DK/DM/DC) 
end 
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for couiU= 1 :leiigtli(ow) 
if meters 

waitbar(coimt/lengtli(ow)); 

end 

z_anal_red=kstat-ow(count)^2 *mstat; 
h_aiial_red=inv(z_anal_red); 

z_exp=k_exp+j *ow(count)*c_expK)w(couiit)^2 *m_exp; 

h_exp=inv(z_exp); 

h_exp=h_exp(aset,aset); 

hacc=h_anal_red(cset_rel,cset_rel); 

lixcc=h_exp(cset__rel,cset_rel); 

dz=mv(inv(iuv(liacc)*(liacc-lixcc)*inv(hacx))diacc); 

ODZ=[ODZ; dz]; 

temp=[temp; eye(cset_size) -ow(count)^2*eye(cset_size) +j*ow'(count)*eye(cset_size)]; 
end 

if meters 

close(waitbar_handle) 

end 

DKDMDC 1 =inv(temp'*inv(R)*temp)*teinp'*inv(R)*ODZ; 

%DKDMDCl=tanp\ODZ ; 

ODKl=DKDMDCl(l:cset_size,:) ; 

ODMl=DKDMDCl(cset_size+l:2*cset_size,:) ; 

ODC 1 =DKDMDC 1(2* cset_sizeH 1 : 3 *cset_size, : ) ; 

clear hacc lixcc h_exp z_exp z_anal_red h_anal_red dz 
clear temp ODZ 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%% 
load stat; 
odcorrha=[]; 

kstat(cset_rel,cset_rel)=kstat(cset_rel,cset__rel>K)DK 1 ; 
mstat(cset_rel,cset_rel)=mstat(cset_rel,cset_rel)+ODM 1 ;; 
cstat(cset__rel,cset_rel)=cstat(cset_rel,cset__rel)+ODC 1 ; 

if meters 

waitbar_haiidle=waitbar(0, ‘Computing Corrected FRF); 
else 

disp('Getting DK/DM/DC') 
end 

for i=l:length(owref) 
if meters 

vvaitbar(i/lengtli(owref)); 

atd 

tempza=kstat+j *owref(i)*cstat-owref(i)^2 *mstat; 
tempha=iiiv(tempza); 
odcorrha=[odcorrha tempha(:)]; 
end 

if meters 

close(vvaitbar_liandle) 

end 

if loopindex = 1 
odcorrha 1 =odcorrha; 
save odl odcorrha 1 
clear odcorrha 1 
elseif loopindex =- 2 
odcorrha2 =odcorrha; 
save od2 odcorrha2 
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clear odcorrha2 
elseif ioopiiidex == 3 
odcorrha3=odcoiTha; 
save od3 odcorrha3 
clear odcoirha3 
elseif loopindex == 4 
odcorrha4=odcorrha; 
save od4 odcorrha4 
clear odcorrha4 
elseif Ioopiiidex == 5 
odcorrha5=odcoiTha; 
save odS odcorrhaS 
clear odcorrhaS 
end 

clear tempha tempza kcorrected mcorrected ccorrected 
clear ODMl ODCl DKDMDCl ODZ temp kstat mstat 

load uiicorrha 
loadH^T^EXP 
if complete 

plotwliichcoord=[cset(niid_index)]; 

else 

plotv^liichcoord=[cset(mid__index)]; 

end 

for plotindex=l :lengtli(plotw^liichcoord) 
which_coord=plotwliichcoord(plotmdex); 
h=figure(h+l); 
fignuni=fignum+l ; 
subplot(2,l,l); 
plot(owret7(2 *pi),. . , 

logl0(abs(H_T_EXP(ndx3d([aset_si2e aset_size lengtli(owref)],.,. 
which_coord,wiiich_coord,':')))),V- . . 
owref7(2*pi),... 

Iogl0(abs(odcorrha(ndx3d([aset_size aset_size lengtli(owre0],,,, 

wliich_coord,wliich_coord,V)))),*g~ « • 
owref/(2*pi),,.. 

Iogl0(abs(uncorrha(ndx3d([aset_size aset_size length(ovvref)],... 
which^coordjWliich^coord,':’)))),^.*) 
ylabel( [^(',int2str(aset( wliich_coord)), ... 

V,int2str(aset(which_coord)),') in/lbf(logl0 of)’]) 
hold on 
v=axis; 

%xlabel('Omega (Hz)') 
if titles 
if use_antires 

title([int2str(odivisions),' PT MATRIX SOLTNS WITH ANTIRESONANCES WEIGHTED LEFf TO RIGHT]) 
else 

tiUe([mt2str(odivisions); PT MATRIX SOLTNS WEIGHTED LEFT TO RIGHT]) 
end 
end 

plot(resfreq( 1 :loopindex)/(2*pi),.. . 
ones(l,loopindex)*(v(3)+abs(v(4)-v(3))*.98),'k+') 

% if use_antires 

% %plot(diag(antiresfreq( 1 :loopindex,l : loopindex ))/(2*pi),ones( 1 ,loopinde%x)*(v(3>i-abs(v(4)-v(3))*.98),'k*') 

% end 
grid on 

lili=legendCExperimental FRFVCorrected MAT Anal FRF,'Uncorrected Anal FRF,Tncluded Modes’); 

axes(lili) 

hold off 
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%if pswitch==y 
% prtfigl(fignum) 

% delete(h) 

% % else 

% delete(h) 

% end 

subplot(2,l,2) 

Qietop=find(abs(owref-ow(lengtli(ovv))* 1 .5 ♦ones( 1 ,iengtli(owreO))==niin(abs(o\vref- 
ovv'(lengtli(ovv))* 1 .5*ones( 1 ,Iengtli(owreO)))); 

plot(owref( 1 :tlietop)/(2*pi),. . . 

iogl0(abs(H_T_EXP(ndx3d([aset_size aset_size lengtli(ovvreO],... 
which_coord,wliicli_coord, 1 :tlietop)))),’r-',.- 
ovvref( 1 :thetopy(2*pi),. . . 

loglO(abs(odcorTha(ndx3d([aset_size aset^size iengtli(ovvTeO],.- 
wliich_coord,which_coord, 1 :thetop)))),’g— 
owref( i :thetop)/(2*pi),... 

ioglO(abs(imcoiTha(ndx3d([aset_size aset_size length(owreJO],... 
wliich_coord,wliich_coord, 1 :tlietop)))),V.’) 

ylabel(['H(’,mt2str(aset(wliich_coord)),V,int2str(aset(wliich_coord)),’) in/lbf (log 10 oO']) 

hold on 

v=axis; 

xlabel(’Omega (Hz)') 
if titles 
if use^antires 

title([mt2sti(odivisions),' PT MATRIX SOLTNS WITH ANTIRESONANCES WEIGHTED LEFT TO RIGHT]) 
else 

tiUe([int2str(odivisions),’ PT MATRIX SOLTNS WEIGHTED LEFT TO RIGHT]) 
end 
end 

plot(resfreq(l:loopindex)/(2*pi),ones(l,loopindex)*(v(3)+abs(v(4)-v(3))*.98),'k+') 

% if use_antires 

% %plot(antiresfreq(which_coord,l:loopindex)/(2*pi),ones(l,loopindex)*(v(%3)+abs(v(4)-v(3))*.98),'k*') 

% end 

grid on 

%lih=legendCExperimental FRF','Corrected MAT Anal FRF,’Uncorrected Anal FRF’j’Included Modes'); 

%axes(hli) 
hold off 

if loopindex == 2 
print -dcdjcolor 
print -dinfile fig6_l 
elseif loopindex == 3 
if odivisions == 1 
print -dcdjcolor 
if complete 
print -dmflle fig6_9 
else 

print -dmflle fig6_6 
end 
end 

elseif loopindex == 4 
print -dcdjcolor 
print -dmflle fig6_2 
end 

% if pswitch==y 

% prtfig 1 ( fignum) 

% delete(li) 

% else 
% delete(h) 

% end 
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380 end ' 
delete(h) 
clear H_T_EXP 
clear odcorrha 
clear uncoirha 
385 end 
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%%%%CHAP6_.2.M%%%%%%%%%%%%%%%%%%%%%%%%%%%%0/00/00/^,0/^0/^0^0^0^0^0/^0^0^^^^^^ 

%%%%%%%%% 

%%% Decompose DZ into DM, DK, & DC %%%%%%%%%%%%%%%%%%%%%%%%%%o/o 0 /o%%%o/o 0 /o%o/o 0 /o% 
%%% using matrix fonnulation %%%%%%%%%%%%%%%%%%%%%%%%%%%%%o/o%%%o/o%o/o% 

%%%%%%%%%%%%%%%%%%%%0/oO/oO/oO/oO/^0/^0/^0^0^0^0^0^0^0/^0/^0/^0^0^0/^0^0^0/^0^0^ 

%%%%%%%%%%%% 

clc 

clear 

closeall 

loadINT 

whitebgCwhite') 

close 

h=0; 

fignum=l; 

iise_antires=l 

titles=0 

mid_index=romid(lengtli(cset)/2); 

clear L Z_ANAL c_exp conn h_anal_red 
clear kexto m_exp mexto tenipj_diags 

clear true_damping true_mass tnie_stifiiiess z_anal z_anal_red z_exp 
clear L Z_ANAL c__exp comi h_anal_red 
clear kexta kexto mextO temp_l_diags 

%%%%%%%%%%%%%%%%%%%FORCETHECSET%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%SPATIALLY INCOMPLETE%%%%%%%%%%%%%%%%%% 
if complete 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%SPATIALLYCOMPLETE%%%%%%%%%%%%%%%%% 
cset=[4 5 6 7 8 9 10 1 1]%%%%%%%%%%%%%%%%%%%%% 
cset_rel=cset%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
cset_size=length(cset)%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
else 

%%%%%%%%%%%%SPATIALLY INCOMPLETE%%%%%%%%%%%%%%%%%% 
cset=[l 3 5 7 9 11 13 15]%%%%%%%%%%%%%%%%%%%%%% 
cset_,rel=[l 2 3 4 5 6 7 8]%%%%%%%%%%%%%%%%%%%%% 
cset_size=lengtli(cset)%%%%%%%%%%%%%%%%%%%%%%%%% 
end 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%SPATIALLY COMPLETE%%%%%%%%%%%%%%%%% 

%cset=[7 8 9 10 11 12 13 14]%%%%%%%%%%%%%%%%%%%%% 
%cset_rel=cset%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%cset_size=lengUi(cset)%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

load exp 
odivisioiis=3 ; 
load odconf 
nunmiodes=9 ; 
odlengtli=2*pi ; 
firsttiine=l; 

dw=odlengtli/(odi visioiis+ 1 ); 
dovv=(freqtop-freqbottom)/(fmeiiess-l ); 
owref= freqbottom:dow:freqtop; 

if firsttiine 
if meters 



127 



450 

455 

460 

465 

470 

475 

480 

485 

490 

495 

500 

505 



CHAP6.M 



waitbarJiandle=vvaitbar(0, 'Computing Experimental FRF'); 
else 

disp('Getting experimental FRF’) 
end 



H_T_EXP=[]; 
for count=l:lengtli(owref) 
if meters 

waitbar(comit/lengtli(owref)); 

end 

*owref(coimt)*c_exp-owref(comit)^2 *m_exp; 
h__exp=inv(z_exp); 
h_exp=h_exp(aset,aset); 

% skyindex=l; 

% for indexl=l :aset_size 
% for index2=index 1 :aset_size 

% red_liolder(skyindex)=h_exp(index 1 ,index2); 

% skyindex=skyindex+ 1 ; 

% end 

% end 

% H^T_EXP=[H^T_EXP redjiolder’]; 

H^T_EXP=[H_T_EXP h_exp(:)]; 
end 



if meters 

close(waitbar_handle) 

end 



save H_T_EXP H_T_EXP 
clear H_T_EXP 
end 

%load H_ANAL^RED; 

%uncorrha=H_ANAL_RED; 

%clear H^ANAL^RED 
if firsttime 
micorrha=[]; 
if meters 

waitbar_handle=waitbar(0,'Computing Uncorrected FRF'); 
else 

dispCGetling Uncorrected FRF') 
end 

load stat 

for i= 1 :lengtli(owre0 
if meters 

waitbar(i/lengtli(owreO); 

end 

tempza=kstat+j *owref(i)*cstat-owref(i)^2 *mstat; 
tempha=inv(tempza); 
micorrha=[uncorrha tempha(:)]; 
end 

if meters 

close(\vaitbar_handle) 

end 

save uncorrha micorrha 
clear uncorrha 
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end 

load H^T_EXP 
for index=l:nummodes 
for coord=l:aset_size 
if index == 1 

antiresffeqs=find(owref>0 & owref<omegax( index)); 
else 

aiitires£feqs=fmd(owref^omegax(index-l) & owref<omegax( index)); 
end 

antires=niin(logl0(abs(H_T_EXP(ndx3d([l aset_size*(aset_size+l )/2 length(o\vreO],... 
l,skyred(coord,coord),antiresffeqs))))); 

wlire=fmd(logl0(abs(H_T_EXP(iidx3d([l aset_size*(aset_size+l)/2 lengtli(owref)],... 

l,skyred(coor(i,coord),antiresfreqs))))==antires); 
temp=owref(antiresfreqs) ; 
antiresfreq(coord,index)=temp(whre); 
end 
end 

for index- l :nummodes 

resffeqs=fmd(owref>oinegax(index)-dow & owref<omegax(iiidex>fdo\v); 
res=ma\(logl0(abs(H_T_EXP(ndx3d([aset_size aset size length(ovvref)]v- 
cset_rel( 1 ),cset_rel( 1 ),resfreqs))))); 

wlire=fmd(logl 0(abs(H_T_EXP(ndx3d([aset_size aset_size length(owref)],... 

cset_rel( 1 ),cset_rel(i ),reslfeqs))))==res); 
temp=o\vref(resfreqs) ; 
resfreq(index)=temp(wlire); 
end 

clear H_T_EXP 
%clear omegax 
%omegax=resfreq; 

for loopindex=startloop:skiploop:endloop 
subdivisions=9; % should be odd for sinipson’s rule 
intlengtli=2*pi; % 1 Hz bandwidtli 
dW=intlengtli/subdivisions; % sanipling freq =.25 HZ 
W=[]; 

W1=G; 

\veiglit=[]; 

lowwei^owref( 1 )-.5 *odlengtli; 
uppper=owref( 1 .5*odlength; 

W=[W lo\vwerHiW:dW:uppper-dW]; 

W 1 = [W 1 lomverf d W :d W;uppper-dW]; 
weiglit=[weiglit ones(size(lowwert-dW:dW:uppper-dW))]; 

for index= 1 iloopindex 
lo\v\ver=omegax( index)-. 5 * intleng tli; 
uppper=omegax(index)+. 5 * intlengtli; 

W=[W low\verH“dW:dW:uppper-dW]; 

W1=[W1 lomverHiW:dW:uppper-dW]; 
weight=[weiglit ones(size(lowwer+dW:dW:uppper-dW))]; 
lo\ver=antiresff eq(index)-. 5 * intlengtli; 
upper=antiresffeq(index>i-.5*intlength; 

W=[W lowwer+dW:dW:uppper-dW]; 

W1=[W1 lowwerfdW:dW:uppper-d^; 
weight=[\veight oues(size(lowweri-dW:dW :uppper-dW))]; 
end 

%%coinpute integrals%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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DZ.R=0; 

DZJ=[]; 
if meters 

waitbar_haiidle=waitbar(0,'Computiiig DZ*); 
else 

dispCComputing DZ’); 
end 

load stat 

for count=l:lengtli(W) 
if meters 

\vaitbar(count/length(W)); 

end 

z_anal_red=kstat/(j *W(count))^-cstat/(j * W(count))+-j * W(comit)*mstat; 
h_anal_red=iiiv(z_anal_red); 
hacc=li__anal_red(cset_rel,cset_rel); 

z_exp=k__exp/(j*W(coimt)>i-c_exp/(j*W(count))+-j*W(count)*m_exp; 

h_exp=inv(z_exp); 

h_exp=h__exp(aset,aset); 

lixcc=h_exp{cset_rel,cset_rel); 

dz=inv(inv(inv(hacc)*(hacc-hxcc)*inv(hacc))-hacc); 

dz_r=real(dz); 

dz_i=imag(dz); 

DZ__R=PZ_R dz_r(:)]; 

DZJ=pZJdzJ(:)]; 

end 

if meters 

close(waitbarJiandle) 

end 

W1=1.AV1; 

[INTK, INTM, INTC]=mtsub(DZ_I,DZ_R,W,Wl,weiglU,cset_size); 
iiitcorrha=[]; 

clear hacc lixcc h_exp z_exp z_anal_red h_anal_red dz 
clear temp ODZ 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%% 
load stat; 
odcorrha=[]; 

kstat(csetJ■el,cset_rel)=kstat(cset_rel,cset_rel)^-INTK; 

instat(cset_rel,cset_rel)=instat(cset_rel,cset_rel)^-INTM; 

cstat(cset_rel,cset_rel)=cstat(cset_rel,cset_rel)+-INTC; 

if meters 

waitbar_handle=\vaitbar(0,’Computing Corrected FRF’); 
else 

disp(’Getting DK/DM/DC’) 
end 

for i= 1 :lengtli(owref) 
if meters 

waitbar(i/leiigtli(o\vref)); 

end 

tempza=kstat+j * o>^ref( i ) * cstat-owref( i)^^2 *mstat; 

tempha=inv(tempza); 

odcorrha=[odcorrha tempha(:)]; 
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end 

if meters 

close(waitbar_haiidle) 

end 



if Icwpindex == 1 
lodcorrhal =odcorrha; 
save lodl lodcorrhal 
clear lodconrhal 
elseif loopindex = 2 
lodcoirha2=odcorrha; 
save lod2 lodcorrhal 
clear lodcorrhal 
elseif loopindex == 3 
lodcorrhaS =odcorrha; 
save lod3 lodcoirhaS 
clear lodcorrlia3 
elseif loopindex = 4 
lodcorrha4=odcoirha; 
save lod4 lodcorrha4 
clear lodcoixha4 
elseif loopindex == 5 
lodcoirha5=odcorrha; 
save lodS lodcorrhaS 
clear lodcorrhaS 
end 



clear tempha tempza kcorrected mcorrected ccorrected 
clear INTM INTC INTK temp kstat mstat cstat 
load micoirha 
load H_T_EXP 

if complete 

plotwliichcoord=[cset(nnd_mdex)]; 

else 

plotwliichcoord=[cset(nnd_index)]; 

end 

for plotindex=l:length(plotwliichcoord) 
wliich_coord=plotwhichcoord(plotindex); 
h=rigure(h+l); 
subplot(2,l,l) 
plot(owref/(2 *pi),. . . 

logl0(abs(H_T_EXP(ndx3d([aset_size aset_size lenglli(owref)],.., 
wliich_coord,wliich_coord,’:’) ))),’r-’, . . . 
o\vref/(2^pi),.., 

Iogl0(abs(odcorrha(ndx3d([aset_size aset_size lengtli(ovvref)],... 
wliich_coord,wliich_coord,V)))),’g— . . 
o>vref7(2*pi),... 

log 1 0(abs(uncorrha(ndx3d( [aset_size aset_size lengtli(o wref)],- . . 
wliich_coord,wlnch_coord,’;')))),1>.') 
ylabel(['H(',intlstr(aset( which^coord)),. . . 

'/,intlstr(aset( wliich_coord)),') in/lbf (log 1 0 of)'] ) 
hold on 
v=axis; 

%xlabel(’Omega (Hz)') 
if titles 

if use antires == 'on ' 
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title([int2str(cxlivisions),' PT INTEGRAL SOLTNS WITH ANTIRESONANCES WEIGHTED LEFT TO RIGHT]) 
else 

Utle([int2str(odivisions),' PT INTEGRAL SOLTNS WEIGHTED LEFT TO RIGHT]) 
end 
end 

plot(oniega\( I :lcx)pindex)/(2*pi),... 
ones(lJoopmdex)*(v(3>fabs(v(4)-v(3))*.98),'k+') 

% if use_antires 

% %plot(aiitiresfreq(wliich_coord,LToopindex)/(2*pi),ones(I,loopindex)*(v(%3)+abs(v(4)-v(3))*.98),’k*') 

% end 
grid on 

hh=legend(Experimental FRF’,’Corrected INT Anal FRF’,TJncorrected Anal FRF,'Included Modes'); 

axes(hli) 

hold off 

%if pswitch=='y' 

% prtfig 1 (fignuin) 

% delete(h) 

%else 

% delete(h) 

%end 

subplot(2,l,2) 

thetop=fmd(abs(owref-W(length(W))* 1 .3*ones( 1 ,length(owreO))=min(abs(owref- 
W(lengUi(W))* 1 .3 ♦onesC Uengdi(owref))))); 

% h=figure(h+l ); 

% fignum=fignum+l ; 
plot(owref(l :thetop)/(2*pi),... 

logl0(abs(H_T_E^(ndx3d([aset_size aset_size length(owre0],... 
wliich_cx)ord,wliich_coord, 1 :tlietop)))),'r-',... 
owref( 1 :tlietop)/(2*pi),. . . 

Iogl0(abs(odcorrha(ndx3d([aset_size aset_size IengtIi(o\\Te0],... 
wliich_coord,wliich_coord, 1 :tiietop)))),’g-',... 
owref( 1 :thetop)/(2*pi),. . . 

loglO(ab$(iincorTha(ndx3d([aset_size aset_size IengUi(o\vre0],-.- 
wliich_coord,wliich_coord, 1 itiietop)))),^.') 

ylabel(pH(',int2str(aset(wliich_coord)),',',int2str(aset(\vIiich_coord)),') in/lbf (loglO oQ’]) 

hold on 

v=axis; 

xlabel('Omega (Hz)') 
if titles 

if use_antires == 'on ' 

Utle([mt2str(odivisions),’ PT INTEGRAL SOLTNS WITH ANTIRESONANCES WEIGHTED LEFT TO RIGHT]) 
else 

title([int2str(odivisions),' PT INTEGRAL SOLTNS WEIGHTED LEFT TO RIGHT]) 
end 
end 

plot(omegax( 1 :loopindex)/(2*pi),ones( 1 ,loopindex)*(v(3)+abs( v(4)-v(3))*.98),'k+') 

% if use_antires 

% %plot(antiresfreq(which_coord, 1 :loopindex)/(2*pi),ones( l,Ioopindex)*(v(%3)+abs(v(4)-v(3))*.98),'k*') 

% end 
grid on 

%lili=legend(Experiineiital FRF','Corrected ENT Anal %RF','Uncorrected Anal FRF','hicluded Modes'); 

%axes(lili) 



if loopindex == 2 
print -dcdj color 
print -dmfile fig6_3 
h=figure(h+l); 
hold off 
load od2 

pIot(ovvref( 1 :tlietop)/(2 *pi),... 
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logl0(abs(H_T_EXP(udx3d([aset_size aset__size lengUi(o\vref)],.- 
\viiich_coord, wliich_coord, 1 :tlietop)))),'r-’,. . . 
owref( 1 :tlietop)/(2*pi),... 

Iogl0(abs(odcorrha(ndx3d([aset__size aset_size lengtli(o\vreO], ••• 
v\iiicli_coord,\vliich_coord, 1 :tlietop)))),’g~ . . 
owref( 1 :tlietopy(2*pi),. .. 

loglO(abs(odcoiTha2(ndx3d([aset_size aset_size lengtli(owref)],- * 
wliich__coord,vvliich_coord, 1 :tlietop)))),’b-.',. . . 
owref( 1 :tlietopy(2*pi),. . . 

logl0(abs(iincx)iTha(ndx3d([a$et_size a$et_size length(owreO],. - 
which_coord,which_coord, 1 :tlietop)))),’m:’) 

ylabel([H(’,int2str(aset(wluch_coord)),’,\mt2str(aset(wliich_coord)),') in/lbf (loglO of)’]) 

hold on 

v=axis; 

xlabel(’Omega (Hz)’) 
if titles 

if use_antires == 'on ’ 

tiUe([int2str(odivisions); PT INTEGRAL SOLTNS WITH ANTIRESONANCES WEIGHTED LEFT TO RIGHT]) 
else 

title([int2str(odivisions),' PT INTEGRAL SOLTNS WElGPfTED LEFT TO RIGHT]) 
end 
end 

plot(omegax(l:loopiiidexy(2*pi),ones(l,loopindex)*(v(3)+abs(v(4>-v(3))*.98),'k+’) 

% if use_antires 

% %plot(antiresffeq(vvliich_coord, 1 :loopmdex)/(2*pi),ones( 1 ,loopindex)*(v(%3)f abs( v(4)-v(3 ))* .98),'k*') 

% end 

grid on 

lili=legend(Experiniental FRF’,'Corrected INT Anal FRF,'Coirected MAT Anal FRF',’Uncorrected Anal FRF,'Includcd 
Modes'); 

axes(lih) 
hold off 
print -dcdjcolor 
print -dinfile fig6_5 
elseif loopindex — 3 
if odivisions == 1 
print -dcdjcolor 
print -dnifile fig6_7 
h=figure(h+l); 
hold off 
loadod3 

plot(owref( 1 :tlietop)/(2 *pi),. . . 

logl0(abs(H__T_EXP(ndx3d([aset_size a$et_size length(ovvre0],... 
which_coord,wliich_coord, 1 :tlietop)))),'r-',., . 
o\vref( 1 :llietopy(2*pi),... 

logl0(abs(odcoiTha(ndx3d([aset_size aset_size lengtli(o\\Te0],. * 

\\'iiich_coord,wliich_coord, 1 :tlietop)))),'g— . . 
owref( 1 :tlietopy(2*pi),.. . 

Iogl0(abs(odcorrha3(ndx3d([aset_size aset_size lengtli(owref)],... 

which_coord,\vliich_coord, 1 itlietop)))),!)-.',. . . 
ov\xef( 1 :tlietop)/(2*pi),. . . 

Iogl0(abs(unconrha(ndx3d([aset_size aset_size lengtli(o\vref)],... 
wliich_coord, which_coord, 1 :tlietop)))),'m:') 

ylabel([’H(’,int2str(aset(wliich_coord)),',',int2str(aset(\vhich_coord)),’) in/lbf (log 10 of)']) 

hold on 

v=axis; 

xlabel('Oniega (Hz)’) 
if titles 

if use__antires -- 'on ' 

tiUe((int2str(odivisions),' PT INTEGRAL SOLTNS WITH ANTIRESONANCES WEIGHTED LEFT TO RIGHT]) 
else 
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tiUe([int2str(odivisions),' PT INTEGRAL SOLTNS WEIGITTED LEFT TO RIGHT]) 
end 
end 

plot(omegax( 1 :loopindex)/(2*pi),ones( 1 ,loopindex)*(v(3)+abs(v(4)-v(3))*.98),’k+') 

% if use_antires 

% %plot(antiresfreq(\vliich_coord, 1 :loopindex)/(2*pi),ones( 1 ,loopindex)*(v(%3)+abs(v(4>-v(3))*.98),'k*') 

% end 

grid on 

lili=legendCExperiniental FRF’,'Corrected INT Anal FRFVUncorrected Anal FRl"',’Included Modes'); 

axes(lih) 

hold off 

print -dcdj color 

print -dmfile fig6_8 

end 

elseif loopindex == 4 
print -dcdj color 
print -dmfile fig6_4 
end 

% if pswitch==’y' 

% prtfigl(fignum) 

% delete(h) 

% else 
% delete(h) 

% end 
end 

delete(h) 
clear H_T_EXP 
clear odcorrha 
clear uncorrha 
end 

loadH_T_EXP 
load odl 
load od2 
loadod3 
load od4 
load uncorrha 
%load lod5 

h=figure(h+l); 

plot(owreI7(2*pi),... 

loglO(abs(H_T_EXP(ndx3d([aset_size aset^size lengtli(owref)]v- 
which_coord, wliich_coord,':')))),'r- , . . . 
o\\Tef7(2*pi),... 

Iogl0(abs(odcorrhal(ndx3d([aset_size aset_size lengtli(owre0]v 
\vliich_coord,\vhich_coord,':')))),’g— . . 
owref7(2*pi),... 

Iogl0(abs(odcorrha2(ndx3d([aset_size aset_size length(ov\Tef)],... 
which_coord,which_coord,':')))),'b-.',. . . 
owTef7(2*pi),... 

Iogl0(abs(odcorrha3(ndx3d([aset_size aset_size length(owTe0]v 
\vliich_coord,wliich_coord,’:')))),'m— . . 
o\vref/(2*pi),... 

Iogl0(abs(odcorrha4(ndx3d([aset_size aset_size length(owref)],... 
wliich_coord,\vliich__coord,':')))),'c-. . 
o\\Tef7(2*pi),... 

Iogl0(abs(uncorrha(ndx3d([aset_size aset_size length(o\\Te0], ••• 
which_coord,wluch_coord,':')))),’k-.') 
yiabel([’H(’,int2str(aset(\vliich_coord)),... 

',’,int2str(asel(wliich_coord)),') in/lbf(logl0 oO']) 
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hold on 
v=axis; 

%xlabel(’Oinega (Hz)') 
if lilies 

if use_antires == 'on ' 

litle([int2slr(odivisions),' PT INTEGRAL SOLTNS WH'H ANTn^SONANCES WEIGITEED LEFT TO RIGHT]) 
else 

title([iiit2slr(odivisions); PT INTEGRAL SOLTNS WEIGHTED LEFT TO RIGITT]) 
end 
end 

% plol(omegax(l :loopindex)/(2*pi),... 

% ones(l,loopiiidex)*(v(3)+abs(v(4)-v(3))*.98),'k+’) 

% if use_anlires 

% %plol(anliresfreq(vvhich_coord,l:loopindex)/(2*pi),ones(I,loopindex)*(v(%3)+abs(v(4>v(3))*.98),'k*') 

% end 
v=axis; 

axis([10 300 v(3) v(4)]), 
grid on 

lili=legend('Experiinenlal FRF’/l mode MAT solulioii','2 mode MAT solulion',’3 mode MAT solulion',’4 mode MAT 
solution'jUncorrected Anal FRF); 
axes(lili) 
hold off 

%if pswilch==y 
% prtfig 1 (fignum) 

% delete(h) 

%else 

% delete(h) 

%end 



prinl -dcdjcolor 
print -dmfile fig6_12 
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%%%%%%%%%% 

% PROGRAM: SETUP.M % 

% INTITLIZES DATA FILE FOR A FINITE ELEMENT ANALYSIS % 

% VARIABLES SAVED TO FILE SETUP.MAT % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%0/o0/o0/o0/o0/^0/^0/^0/^0/^0/^0/^0/^0/^0/^0/^0/^^^^^^ 

%%%%%%%%%% 

clc 

clear 

L=inputCEnter tlie total lengtli of tlie beam: '); 
iiumel=inputCNum of elements: '); 
dofjiode=iiiput('num dof per node: '); 
area=input(’area of beam: ’); 
eeii=inputCEI of beam: '); 

%ee=iiiput('modulus for beam: '); 
pho=input(’mass density for beam: '); 
numdofmof_node*(numel+ 1 ) 
clc 

wlule 1 

conforce=inputCEnter # of concentrated loads: *); 
if ~isnan(conforce) 
break 
end 
end 

for i=l:conforce 

forcepos(i)=iiiput(['Beam position of force ’,nuin2str(i),' : ']); 
forcesiz(i)=input( [Magnitude of force ’,num2str(i),' : ’]); 
end 

force=zeros(nunidof,l ); 
for i=l:conforce 

temppos=ceil(forcepos(i)/(L/numel)); 

force(dof_node *temppos+ 1 )=force(dof_node* temppos+ 1 5 * forcesiz( i ); 

force(dof_node*temppos+3)=force(dof_node*temppos+3>»'.5*forcesiz(i); 

end 
vvliile 1 

lmass=iiiput('Enter # of Lumped masses: '); 
if ~isnan(lmass) 
break 
end 
end 

for i=l :lmass 

masspos(i)=input(['Beam position of Mass *,num2str(i),’ : ']); 
masssiz(i)=iiiput([Magnitude of Mass '4iuni2str(i),' : ’]); 
end 

lumpinass=zeros(numel, 1 ); 
for i=l:lmass 

temppos=round(masspos(i)/(L/numel)); 
lunipmass(temppos)=lumpniass( temppos)+masssiz( i ); 
end 

lumpmass' 
wliile 1 

lspring=mput('Enter U of Lumped Springs: '); 
if-isnan(lspring) 
break 
end 
end 

for i=l:lspring 

springpos(i)=input(['Beam position of Spring ',iium2str(i),' : ']), 
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springsiz(i)=iiiput(['Sprmg constant of Spring ’,nimi2str(i),' : ']); 
end 

liunpspring=zeros(mmiel,l ); 
for i= 1 lispring 

temppos=round(springpos(i)/(L/nuinel)); 

lunipspring(temppos)=luinpspring(temppos)+springsiz(i); 

end 

lumpspring’ 
conn=[l^]; 
for i=2:numel 
comi=(conn;i,i+l]; 
end 
conn 
while 1 

bc=['pimied-pinned ’ 

’clamp-clamp ’ 

’left guided clamp ' 

’riglit guided clamp’ 

’cantilevered ’ 

’free-free ’]; 
clc 

help bctext 

n=input(’Select a boundary condition: ’); 
if((n>0)&(n<7)) 
break 
end 
end 

ifn == 1 
bc=*pp’; 
elseif n == 2 
bc='cc’; 
elseif n = 3 
bc=lc’; 
elseifn==4 
bc='rc’; 
elseif n==5 
bc=’cl' 
else 
bc=’fF; 
end 

clear i 

clear temppos 
clear conforce 
clear forcepos 
clear forcesiz 
clear Imass 
clear masspos 
clear masssiz 
clear Ispring 
clear springpos 
clear springsiz 
save setup.mat 
x=0:L/numel:L; 
for i=l:lengtli(x) 
h(i)=5; 
end 
clg 

hold off 
%axis('ofr) 
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%plot([],[]) 
%hoId on 
plot(xJi,x,h,V) 
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%%%%%%% 

%FORM FE BEAM MODEL WITH ERRORS AT PRESCRIBED LOCATIONS%%%%%%%%%%o/o 
%MASS ERROR LOCATION SPECIFIED VIA POSOFMASSER. VALUE OF ERROR IN%% 

%PERCENT GIVEN BY VALOFMASSERR, ETC.%%%%%%%%%%%%%%%%%%%%%%o/o%%o/oO/oO/oO/oO/oO/^ 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%0/oO/oO/oO/^0/^^0^0/^0/^0/^0^0^0/^0/^0^0^0^0/^% 

%%%%%%% 

clear,clc;clg 

casenanie-Spatially Incomplete A set Stiffness error'; 
if exist('setup.mat') 

load setup %load existing input data file 
else 

setup %create new input data file 

end 

slnic_damping=0.0000000000 1 ; 

Iunipdamp=lumpspring*0 ; 
posolMassen=[3 4] ; 
valueofMasserr=[0.25 0.25] ; 
posofStiffen^(3 4j; 
valueofStifferr=[0.25 0.25] ; 
posofDamperr=[3 4]; 
valueofDaniperr=[0.25 0.25] ; 
for i=l :lengtli(posofMasserr) 
luinpmass(posofMasserr(i))=valueofMasserr(i); 
end 

for i=l:lengtli(posofStifferr) 
lunipspiing(posofStifferr(i))=valueofStifreiT(i); 
end 

for i=l:lengtIi(posofDamperr) 
lumpdanip(posofDampeir(i))=valueofDamperr(i); 
end 

elen=L/nuinel; 
ifbc=='pp’ 
doftoldll=2; 
elseif bc='cc’ 
doftokill=4; 
elseif bc=-fT 
doflokill=0; 
elseif bc=- cl' 
doflokill=2; 
end 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%BUILD THE ELEMENTAL MASS AND STIFFNESS MATRICES%%%%%%% 

%FOR A2DOF/NODE "F.E." STRUCTURES. %%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



%ELEMENT STff FNESS MATRIX%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
ke=[12 6*elen-12 6*eleii; 

6*elen 4*(elen^2) -6*elen 2*(elen^2); 

-12 -6*elen 12 -6*eleii; 

6*elen 2*(eleiY2) -6*elen 4*(elen^2)]; 

ke=(eeii/(elen'^3 )). *ke; 
g=386; 

elemke=eeii/(elen^3 ); 

%ELEMENT MASS MATRICE%%%%%%%%%%%%%%%%%%%%%%%%%% 
me=[156 22*elen 54 -13*elen; 

22*elen 4*(elen'^2) 13*elen -3*elen^2; 
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54 . 13*elen 156 -22*elen; 

-13*elen-3*(elen'^2) -22^elen 4*(elen'^2)]; 

me = (((pho*area)*elen/g)/(420)).*me; 
elemnie=(((pho*area)*elen/g)/(420)); 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%% 

%ASSEMBLE GOBAL STIFFNESS AND MASS MATT^ FOR A 2 DOF F.E. STRUCTURE%%%% 

%BASED ON THE ELEMENTAL MATRIXES KE AND ME. THE STRUCTURE CONSISTS OF%%% 

%NUMEL ELEMENTS WITH ELEMENT CONNECTIVITY GIVEN BY MATRIX CONN%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%% 
goblk=zeros(niimdo0; 
goblm=zeros(numdo0; 
goblc=zeros(numdof); 
for i=l:nuniel 
v=comi(i,l); 

\v=comi(i,2); 

goblk(dof_node*v-l :dof_node*v,dof_node*v-l :dof_node*v)=... 
goblk(dof_node* V- 1 :dof_node * v,dof_node* v- 1 :dof_node* v)+. . . 
ke(l :dof_node, 1 :dof_node); 

goblk(dof_node*v-l :dof_node* v,dof_iiode*w-l :dof_iiode*w)=... 
goblk(dof_node* V- 1 :dof_node* v,dof_iiode* w- 1 :dof_node*\v)f . . . 
ke( 1 :dof_iiode,dof_nodei-l :2*dof__node); 

goblk(dof_node* w-1 :dof jiode* w,dof_node* V- 1 idofjiode* v)=. . . 
goblk(dof_node*w- 1 :dof_node* w,dof_node* v- 1 :dof_node* v)+. . . 
ke(dof_node+ 1 :2 *dof_ii(^e, 1 :dof_node); 

goblk(dof_node* w- 1 :dof_node* w,dof_node*\v- 1 :dof_node* vv)= . . . 
goblk(dof_node*w-l :dof_node* w,dof_node* w-1 :dof_node*w)+. . . 
ke(dof Jlode^- 1 :2 * dof_ii(^e,dof_node^- 1:2* dof_node ); 

gobbn(dof_node* v- 1 :dof_node* v,dof_node* v- 1 :dof_node* v)=. . . 
goblin(dof_node* v- 1 :dof_node* v,dof_node* v- 1 :dof_iiode* v)+.. . 
me( 1 :dof_node, 1 :dof_node); 

goblm(dof_node* v- 1 :dof_node* v,dof_node* w- 1 :dof_node* w)=. . . 
goblm(dof_node* v- 1 :dof_node* v,dof_node* w- 1 :dof_node* w)^.. . 
me( 1 :dof_node,dof_node+l :2*dof_node); 

goblm(dof_node*\v- 1 :dof_node* w,dof_iiode* v- 1 :dof_iiode* v)=. . . 
goblm(dof_node* w-1 :dof_node* w,dof_node* v-1 :dof_node* v)^-. . . 
me(dof_node+ 1 :2 *dof_node, 1 :dof_node); 

goblm(dof_node*w-l :dof_iiode*w,dof_node*w-l :dof_node* w)=.. . 
goblm(dof_node*w-l :dof_node*w,dof_node* w-1 :dof_node*w)+.. . 
me(dof_node»- 1 :2* dof_node,dof_node+ 1:2* dof_node); 
end 

goblc=sqrt(- 1 )*struc_damping. * goblk; 

goblkx=goblk; 

goblcx=goblc; 

goblnLX=goblm; 



for i=l:numel 
v=comi(i,l); 

\v=comi(i,2); 

goblkx(dof_node* V- 1 :dof_node* v,dof_node* v- 1 :dof_iiode* v)=. . . 
goblkx(dof_iiode* V- 1 :dof_node* v,dof_node* v- 1 :dof_node* v)^. . . 

1 umpspring( i). * ke( 1 :dof_node, 1 :dof jiode); 

goblkx(dof_iiode* V- 1 :dof_node* v,dof_node* w- 1 :dof_node* w)=. . . 
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goblkx(dof_node * V- 1 :dof_node* v,dof_node * w- 1 :dof_node*\v)+. . . 
lumpspring(i).*ke( 1 :dof_node,dof_node+ 1 :2*dof_node); 

goblkx(dof_node*w- 1 :dof_node *w,dof_node * v- 1 :dof_node * v)=. . . 
goblkx(dof_node*w-l :dof_node*w,dof_node*v-l :dofjiode*v)+... 
lumpspring(i).*ke(dof_node+ 1 ;2 *dof_node, 1 :dof_node); 

goblkx(dof_node*w-l :dof_node*w,dof_node*w- 1 :dof_node*w)=... 
goblkx(dof_iiode* w- 1 :dof_node*\v,dof_node*w- 1 :dof_node*vv)+... 
Imnpspring(i). *ke(dof_node+ 1 :2 *dof_node,dof_node+ 1 :2*dof_node); 

gobliiLX(dof_node* v- 1 :dof_node * v,dof_node* v- 1 :dof_node * v)=. . . 
goblmx(dof_node* V- 1 :dof_node* v,dof_node* v-1 :dof_node* v)+. . . 
luinpinass(i). *me( 1 :dof_node, 1 :dof_node); 

goblmx(dof_node*v-l :dof_node*v,dof_node*w-l :dof_node*w)=... 
goblnix(dof_node * V- 1 :dof_node * v,dof_node * w- 1 :dof_node* w)+. . . 
luinpmass(i). *me( 1 :dof_node,dof_node+ 1 :2 *dof_node); 

goblmx(dof_node*w-l :dof_node*w,dof_node*v-l :dof_node* v)=... 
gobbiix( dof_node * w- 1 :dof_node *w,dof_uode * v- 1 :dof_node * v)^-. . . 
luinpmass(i). *me(dof_node+ 1 ;2*dof_node, 1 :dof_node); 

gobliiix(dof_node*w-l :dof_node*w,dof_node*w-l :dof_node*\v)=... 
gobliiLX(dof_node*w-l :dof_node*w,dof_node*w- 1 idofjiode* w)+... 
luinpmass(i).*me(dof_node+l :2*dof_node,dof_node-i-l :2*dof_node); 



goblcx(dof_node*v-l:dof_node*v,dof_node*v-l:dof_node*v)= .. 
goblcx(dof_node* V- 1 :dof_node* v,dof_node* v- 1 :dof_node* v)+. . . 
luinpdamp(i)*sqrt(- 1 )*struc_damping. *ke( 1 :dof_node, 1 :dof_node); 

goblcx(dof_node* V- 1 :dof_node* v,dof_node* w- 1 :dof_node*w)=. . . 
goblcx(dof_node* v-1 :dof_node*v,dof_node*w- 1 :dof_node*w)+.. . 
lumpdainp(i)*sqrt(-l )*stnic_daiuping. *ke( 1 :dof_iiode,dof_node+l :2*dof_node); 

goblcx(dof_node* w-1 :dof_node* w,dof_node* v- 1 :dof_ncxle* v)=. . . 
goblcx(dof_node*w-l :dof_node*\v,dof_node* v-1 :dof_node*v>f ... 
lumpdamp(i)*sqrt(-l )*struc_damping. *ke(dof_node+ 1 :2 *dof_node, 1 idofjiode); 

goblcx(dof_node*w-l :dof_node* w,dof_node’*'w- 1 :dof_node*w)=... 
goblcx(dof_node*w-l :dof_node*\v,dof_ncxle*w-l :dof_node*w)+... 

lumpdamp(i)*sqrt(-l)*stnic_dainping.*ke(dof_iiode+l:2*dof_node,dof_node+l:2*dof_node); 

end 



numdof=nmndof-doflokill; 

[gk,gni,gc,k_anal,m_anal,c_anal]=fixbcs(goblk,goblin,goblc,bc); 

[gLx,giiLX,gcx,k_exp,in_exp,c_exp]=fixbcs(goblkx,goblmx,goblcx,bc); 
save beamdata 



141 



5 

10 

15 

20 

25 



FSTATIC.M 



function [kstat,mstat]=fstatic(k4n,oset,aset) 

aset_size=lengtli(aset); 

kaa=k(aset,aset); 

kao=k(aset,oset); 

koo=k(oset,oset); 

koa=kao*; 

clear k; 

k=[kaa,kaoJcoa,koo]; 

maa=m(aset,aset); 
niao=ni(aset,oset); 
moo=m(oset,oset); 
moa=niao’; 
clear m; 

in=[inaa^nao;inoa,moo]; 



t_static=-koo\koa; 
T_static=[eye(aset_size); t_static]; 

kstat=T_static’*k*T_static; 

nistat=T_static'*m*T_static; 

end 
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% 

luiiction [kirs,inir s] =firs_taiii(k,iii,oset,aset) 

% 

% tliis function returns tlie IRS reduced stiffness 
% and mass matrices, given tlie unreduced couteiparts. 

% Care must be taken tliat tlie aset and oset vectors coirespond 
% witli tlie existing arrangement of k and m. 

% k and m are UNPARTITIONED nialrices. 

% 

aset_size=lengtli( aset); 

% 

kaa=k(aset,aset); 
kao=k(aset,oset); 
koo=k(oset,oset); 
koa=kao’; 
dear k; 

k=[koo,koa;kaoJcaa]; 

% 

maa=m(aset,aset); 
niao=ni(aset,oset); 
moo=m(oset,oset); 
moa=niao'; 
dear ni; 

m=[moo,nioa;niao,iiiaa] ; 

% 

t_static=-koo\koa; 

T_static = [t_static; eye(aset_size)]; 

% 

kstat=T_$tatic'*k*T_static; 

mstat=T_static’*m*T_static; 

% 

tirs=t_static+inv(koo)*(moa+moo*t_static)*inv(nistat)*kstat; 

T_irs=[tirs;eye(aset_size)]; 

% 

kirs=T_irs'*k*T_irs; 

mirs=T_irs'*m*T_irs; 

% 

% end fmiction firs_tam 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%% 

%tliis function returns a vector U containing modal%%%%%%%%%%%%%%%%%%%%% 

%frequencies (rad/sec)^2 in ascending order along%%%%%%%%%%%%%%%%%%%%%% 

5 %witli associated mode shapes%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%% 



1 0 function [u 1 ,lambda,index]=freqmode(k 4 n) 

[u,lambdat]=eig(m\k); 

[lambda, index]=sort(diag(lambdat)); 
ll=zeros(length(k),length(k)); 
for i=l:lengtli(k); 

15 ll(i,i)=lamb^i); 
end 

lambda=diag(ll); 
forj = l:length(k) 
ul(: j) = u(:,indexO')); 

20 end 
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function [r,c] = iidx3d(siz,i j,k) 

%NDX3D Index into 3-D matrix packed in a 2-D matrix. 

% ELEM = NDX3D([M N P],I,J,K) returns tlie element position ELEM 
% of tlie (i j,k) elements of a M-by-N-by-P matrix wliich is 

% stored in a (M*N)-by-P matrix. For example, tlie tliree M-by-N 

% matrices A 1 are packed into a 2-D matrix using 

% A = [A1(:)A2(:)A3(:)]; 

% If lengtli(I) is ni, LENGTH( J) is n, and LENGTH(K) is p, 

% tlien ELEM will be an m-by-n-by-p matrix. 

% 

% [R,C] = NDX3D([M N P],I,J,K) returns tlie row and colmnn 

% position of the (ij,k) element as stored in tlie nomial 
% 2-D matrix of size (M^N)-by-P. 

% 

% To specify all tlie elements along one dimension use 
% For instance, NDX3D([3 5 4],2:3,’:’,3:4) returns tlie 
% elements for tlie 2-by-5-by-2 matrix. 

% 

% See also ELEM3D, MESHGRID, SLICE. 

% Clay M. Tliompson 1 1-3-92 

% Copyriglit (c) 1 992 by Hie MatliWorks, Inc. 

% SRevision: 1.7 $ SDate: 1993/09/03 14:36:52 $ 

if isstr(i), i = 1 :siz(l ); end 
if isstr(j), j = l:siz(2); end 
if isstr(k), k = l:siz(3); end 

if isempty(i) | iseiiipty(j) | isempty(k), r = []; c = []; return, end 

if any( (i<0) | (i>siz(l )) ), errorflndex I out of range.’); end 
if any( O'<0) j ^>siz(2)) ), erroifliidex J out of range.’); end 
if any( (k<0) | (k>siz(3)) ), error(1ndex K out of range.’); end 

ifnargout==2 
[jj,ii] = ineshgrid0‘,i); 
r = ii(:) + GjCH)*siz(l); 
c = k(:); 

else 

[jj,ii3ck] = mesligrid(j,i,k); 
r = ii + (jj-1 )^siz( 1 ) + (kk-1 )*prod(siz( 1 :2)); 

end 
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fiuiction [K,M,C]=intsub(Z_I, Z_R, omega, omega 1 , W, set_size) 

waitbar_handle=waitbar(0,*Computiiig Integrals'); 

j=sqrt(-l); 

integ=omega.^2. *abs( W); 

rint=real(integ); 

iint=imag(integ); 

aa=mytrapz(omega 1 ,rint)+mytrapz(omegal ,iint)*j ; 
integ=(l ./omega.'^2).*abs(W); 
rint=real(integ); 
iint=imag(mteg); 

bb=mytrapz(omegal ,rint)f mytrapz( omega 1 

integ=abs(W); 

rint=real(integ); 

iint=imag(iiiteg); 

cc=mytrapz(omegal ,rint)+mytrapz(oniegal 
for i=l :set_size 
waitbar(i/set_size); 
for k=i:set_size 

M(i,k)=( 1 /(aa*bb-cc^2))*(bb*mytrapz(omegal 
omega. *Z_I(ndx3d([set_size set_size leiigtli(omega )],... 
i,k, 1 :leiigtli(omega))). *abs( W))-cc*mytrapz(omega 1 . 

1. /omega. *Z_I(ndx3d([set_size set_size lengUi( omega)]... 
,i,k,l:lengUi(omega))).*abs(W))) ; 

M(k,i)=M(i,k); 

C(i,k)=(l/cc)*mytrapz(omegal ,Z_R(... 
ndx3d([set_si2e set_size lengtli(omega)],... 
i,k, 1 :lengtli(omega))). *abs(W)); 

C(k,i)=C(i,k); 

K(i,k)=( 1 /(aa*bb-cc'^2))*(cc*mytrapz(omega 1 ,. . . 
omega. *Z_I(ndx3d([set_size set_size lengtli( omega)]... 

,i,k,l :leugtli(omega))).*abs(W))-aa*mytrapz(omega I ,... 

1. /omega. *Z_I(iidx3d([set_size set_size length(omega)]... 
,i,k,l:lengtlTi(omega))).*abs(W))) ; 

K(k,i)=K(i4c); 

end 

end 

close(waitbar_handle) 



146 



5 

10 

15 

20 

25 



MYTRAPZ.M 



function z = lrapz(x,y) 

%TRAPZ Trapezoidal numerical integration. 

% Z = TRAPZ(X,Y) computes tlie integral of Y witli respect to X using 
% trapezoidal integration. X and Y must be vectors of tlie same lengtli, 
% or X must be a column vector and Y a matrix witli as many rows as X. 

% TRAPZ computes tlie integral of each column of Y separately. 

% Tlie resulting Z is a scalar or a row vector. 

% 

% Z = TRAPZ(Y) computes tlie trapezoidal integral of Y assuming unit 

% spacmg between tlie data points. To compute tlie integral for 

% spacing different from one, multiply Z by tlie spacing increment. 

% 

% See also SUM, CUMSUM. 

% Clay M. Tliompson, 10/16/90; Cleve Moler, 1/19/92. 

% Copyriglit (c) 1 984-94 by Tlie MatliWorks, Inc. 

% Make sure x and y are coluimi vectors, or y is a matrix. 

% Trapezoid sum computed witli vector-matrix multiply. 

[iii,ii]=si 2 e(x) ; 
z=0 ; 

for index=l:ni 

yy=y( l+(index-l )*n:n+(index- 1 )*n); 

yy=yy(.-% 

xx=x(hidex,:); 

XX = xx(:); 

z =z+difl(xx)' *(yy(T.n-l) + yy(2:n))/2; 
end 
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